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. This manuscript is a review of the main results obtained in a series of papers involving the present 

I authors and their collaborator J.L. Cardy over the last two years. In our work we have developed 

CN ■ and applied a new approach for the computation of the bi-partite entanglement entropy in 

' massive H-l-dimensional quantum field theories. In most of our work we have considered these 

I ^ I theories to be also integrable. Our approach combines two main ingredients: the "replica trick," 

' and form factors for integrable models and more generally for massive quantum field theory. Our 

7—1 ■ basic idea for combining fruitfully these two ingredients is that of the branch-point twist field. 

i By the replica trick we obtained an alternative way of expressing the entanglement entropy 

r~| ' as a function of the correlation functions of branch-point twist fields. On the other hand, a 

"V* ' generalisation of the form factor program has allowed us to study, and in integrable cases to 
obtain exact expressions for, form factors of such twist fields. By the usual decomposition of 

r~| I correlation functions in an infinite series involving form factors, we obtained exact results for 
the infrared behaviours of the bi-partite entanglement entropy, and studied both its infrared 

7—1 . and ultraviolet behaviours for different kinds of models: with and without boundaries and 

\ backscattering, at and out of integrability. 
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1 Introduction 

Entanglement is a fundamental property of quantum systems. Its most striking consequence is 
the fact that performing a local measurement may affect the outcome of local measurements 
far away; this is one of the main differences between quantum and classical systems. Given 
that our understanding of the physical world is largely based on everyday experience and that 
entanglement seems to contradict such experiences, it is not surprising that the existence of 
quantum entanglement has been a source of controversy and heated scientific debate for some 
time (see e.g. |l] for the famous EPR-paradox). Today entanglement is a well established and 
measurable phenomenon, whose reality was experimentally confirmed in the early eighties by the 
famous experiments of Alain Aspect and collaborators [21 [3], using pairs of maximally entangled 
photons. In the last decades, many applications of entanglement have developed into successful 
fields of research such as quantum computation and quantum cryptography. Entanglement lies 
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also at the heart of other interesting phenomena such as quantum teleportation (see [1| for a 
review) . 

As a consequence of the prominent role of entanglement in quantum physics, there has been 
great interest in developing efficient (theoretical) measures of entanglement. For example, a 
quantity of current interest in quantum models with many local degrees of freedom is the bi- 
partite entanglement entropy [5], which we will consider in this review. In its most general 
understanding, it is a measure of the amount of quantum entanglement, in some pure quantum 
state, between the degrees of freedom associated to two sets of independent observables whose 
union is complete on the Hilbert space. In the cases considered in this review, the quantum 
state will mostly be the ground state |gs) of some extended 1+1-dimensional (1 space + 1 time 
dimension) quantum model, and the two sets of observables correspond to the local observables 
in two connected regions, say A and its complement, A (we will also briefly discuss the general 
technique in the case of excited states). Other measures of entanglement exist, see e.g. OEUHllll 
llOj . which occur in the context of quantum computing, for instance. Measures of entanglement 
are important at a theoretical level, as they give a good description of the quantum nature of a 
ground state, perhaps more so than correlation functions. 

General aspects of the entanglement entropy in extended quantum systems will be discussed 
in [6]. For our present purposes, prominent examples of extended one-dimensional quantum 
systems are quantum spin chains, which model physical systems consisting of infinitely long 
one-dimensional arrays of equidistant atoms characterized by their spin. Their entanglement 
has been extensively studied in the literature [TUfT^fTHlfnifT^fTCIfTTlfTFlfTU]. 

a) 

V 

A : L sites 

!■ 
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Figure 1: The entangled region A of length L of a quantum spin chain, a) the bulk case, b) the 
boundary case. 



b) 



In order to provide a formal definition of the entanglement entropy, let us consider the 
Hilbert space of a quantum model, such as the chain above, as a tensor product of local Hilbert 
spaces associated to its sites. This can be written as a tensor product of the two Hilbert spaces 
associated to the regions A and A: 

n = A^A. (1.1) 

Then the entanglement entropy is the von Neumann entropy of the reduced density matrix 
associated to A: 

Sa = -T^APA^ogpA , Pa = Tr^|gs>(gs|- (1-2) 

We will be interested in analysing the entanglement entropy in the scaling limit of infinite- 
length quantum chains. The scaling limit gives the universal part of the quantum chain behaviour 
near a quantum critical point, which is described by a model of 1+1-dimensional quantum field 
theory (QFT) (which we will assume throughout to possess Poincare invariance). The scaling 



2 



limit is obtained by approaching the critical point while letting the length L of the region A 
go to infinity in a fixed proportion with the correlation length ^ (these lengths are measured in 
number of lattice sites). In this limit, the entanglement entropy is in fact divergent. The way 
the entanglement entropy diverges was understood in |20 ^ [211 [22] . It is controlled by the central 
charge c corresponding to the critical point that is being approached. In general, for every point 
that separates a connected component of A from a connected component of its complement 
A (a boundary point of A), there is a term c/6 log^. We will look at two cases: where the 
quantum chain is infinite in both directions and A is a connected segment, with two boundary 
points (the bulk case), and where the quantum chain is infinite only in one direction, and A is 
a connected segment starting at the boundary of the chain, with only one true boundary point 
(the boundary case) - see Fig. [H In these two cases, the divergent part of the entanglement 
entropy is respectively: 

5r'~|ioge + o(i), 5— ^^^£ioge + o(i). (1.3) 

These divergent terms are not universal (as they depend on the correlation length); the universal 
terms of the entanglement entropy are hidden in the 0(1) part. These universal terms, that 
depend on the proportion to ^ with which A is sent to infinity, are in general described by QFT 
(conformal, massive, etc.). The analysis of the universal terms using CFT techniques in various 
situations is reviewed in [23]. In the present review, we will explain how to use massive QFT 
techniques in order to obtain information about the universal terms. We will now overview the 
main ideas and results of our works on this subject. 



1.1 Overview of ideas and results 

It is known since some time |20^ [211 [22l [2l] that the bi-partite entanglement entropy in the 
scaling limit can be re-written in terms of more geometric quantities, using a method known as 
the "replica trick" . The essence of the method is to "replace" the original QFT model by a new 
model consisting of n copies (replicas) of the original one, in order to use the formula 

Sa = - lim ^Ttj,p% (1.4) 
n^i an 

The trace in this formula is reproduced by the condition that these copies be connected cyclicly 
through a finite cut on the region A. Then, this trace is the partition function Z„(xi,X2) of 
the original (euclidean) QFT model on a Riemann surface Mn,xx,x2 with two branch points, at 
the points xi and X2 in M^, and n sheets cyclicly connected (we will provide more explanations 
about this in the next section). Fig. [2] shows a representation of such a Riemann surface. The 
positions of the branch points correspond to the end-points of the region A in the scaling limit. 
This gives: 

5^(1.1 -x,|) = -lim-fMpi^. (1.5) 

n^i an Z{ 

Here, |xi — 3;2| is the euclidean distance between xi and X2- This formula holds both for our 
bulk and boundary cases. In the boundary case, one of the branch points (say xi) is on the 
boundary of the model, xi = (0,0). The positions xi and X2 are dimensionful positions in the 
QFT model. They are naturally at zero imaginary time xi = (xi, 0), X2 = (x2, 0) (but this is not 
crucial because the euclidean QFT has rotation invariance), and their x-coordinates are related 
to the ratio between the dimensionless region length L and correlation length ^ of the quantum 
chain by |xi — X2I = L/{mS,), where m is the QFT mass scale associated to ^ (whose only role 
here is to provide a dimension). 
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Figure 2: A representation of the Riemann surface M 



n,xi,X2 • 



Naturally, this expression implies that we must analytically continue the quantity 2:2) 
from n S N, where it is naturally associated to Riemann surfaces, to n G [l,oo). The object 
Tr^p^ certainly has a well-defined meaning for any n such that Re(n) > 0. Indeed, is 
hermitian (and has non-negative eigenvalues summing to 1), so that Tr_Ap^ is the sum of the n^^ 
powers of its eigenvalues (with multiplicities). Note that this is an analytic continuation from 
positive integers n to complex n that satisfies the requirements of Carlson's theorem [25], hence 
the unique one that does. The scaling limit of this object is what defines the proper analytic 
continuation of Z„(xi,X2). Finding the correct analytic continuation has been one of the major 
challenges encountered in our work. It is natural to assume, as it has been done before [20] 
and discussed in [26] in the present context, that the two branch points just become conical 
singularities with angle 27rn, the rest of the space being flat. 

As will be described later in more detail, in [27] we showed that there is a way of asso- 
ciating the branch points at xi and X2 to local QFT fields: through branch-point twist fields 
T(xi),T(x2)- These twist fields are defined only in the replica model (not in the original model), 
and are associated to certain elements of the extra permutation symmetry present in the replica 
model. In terms of these fields we showed that, in the bulk case. 



Here (0| • • • |0) denote correlation functions in the n-copy model; the state |0) is the vacuum 
state of the latter. The branch-point twist fields have the CFT normalisation (which we will 
discuss later). The constant Zn, with Zi = 1, is an n-dependent non-universal constant, e is a 
short-distance cut-off which is scaled in such a way that dZyij dn — at n — 1, and, finally, d^i is 
the scaling dimension of the counter parts of the fields T, T in the underlying n-copy conformal 
field theory. 



which can be obtained by CFT arguments [241 [27] and where c is the central charge. The short- 
distance cut-off is related to the correlation length via e = a/(m^) for some dimensionless finite 
non-universal number a. This short distance cut-off takes care of the infinite contributions to 
the partition functions Zn{xi,X2), as compared to the n**^ power of Zi, around the points xi 
and X2 where a branch point lies. There is one contribution of ej^ for each point, corresponding 
to one contribution of c/61og^ to the entanglement entropy. 




^„e2<i„^0|r(xi)f(x2)|0). 



(1.6) 




(1.7) 
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In our most recent work |28j we have generalized this understanding to the boundary case, 
where now the region A extends between the origin xi = 0, where the boundary of the model is 
located, and the x-coordinate X2. With similar arguments we have showed that 

^I^^=Z^e'^-{0\T{x2)\B). (1.8) 

The state \B) is a boundary state, which depends on the particular model and boundary condi- 
tion under consideration. In the context of integrable models, it was introduced in the seminal 
work of Ghoshal and Zamolodchikov [29J. Here, we take it with the normalisation {0\B) = 1. 
The factors Zn and e are defined through similar conditions as in the bulk case. 

In terms of the variables r = |xi — X2I, e, m, and using the CFT normalisation of the branch- 
point twist fields as well as large-distance factorisation of correlation functions, we can re-write 
the logarithmic divergence formula (jl.3p in a more precise fashion. The IR (large-mr) and UV 
(small-mr) leading behaviours of the entropy in the bulk are 

- log(r/e) -I- 0(1) e < r < 

(1 9) 

-- log(em) + [7"^°^^' + 0{{rmy°°) e < < r 
3 

The term [/™°dei jg model-dependent constant, which we computed exactly for the Ising model 

Rising ^ _o.i3i984... (1.10) 

in [57] using QFT methods, reproducing results of ^3] obtained on the lattice. In p.9|) . the 
short-distance cut-off e is of course non-universal, and it is in general hard to evaluate its exact 
relation to the correlation length. But the UV leading behaviour provides an unambiguous QFT 
definition for it, independently of what the correlation length of a particular quantum chain may 
be. Once this definition is taken, the terms that are added to — |log(e) are universal ~ these 
terms constitute the universal part of the entanglement entropy. For instance, the constant 
jjmodei jg ii^deed a universal QFT quantity. Note that the IR leading behaviour is interpreted 
as the saturation of the entanglement entropy at large distances. 

In the boundary case, our definitions of the boundary state further imply that the UV and 
IR leading behaviours of the entanglement entropy are 

-log(2r/e) + V(k) + 0(1) e < r < m"^ 

6 



^boundary 

~ 6 ' 2 



p [/model (1-11) 

-log(em)H hO((rm)~°°) e < m"^ < r 



Again, e is non-universal and hard to calculate. But since we know that the constant JJ^"'^'^^ is 
universal, the IR behaviour in (jl.lip gives an unambiguous definition of £. Then, the terms that 
are added to — ^ log(e) are the universal part of the entanglement entropy, true QFT quantities. 
In particular, the UV behaviour in (II. lip provides a universal definition for V{k). Note that the 
leading asymptotic term at large distance [/m°deiy'2 in the boundary case is just a choice. Once 
this choice is made, the constant V{k) is universally fixed. In particular, our short-distance cut- 
off e here is in general different in the bulk and boundary cases - it is related to the correlation 
length ^ in different ways. Here, k is a parameter that represents the boundary condition. 

A major focus of our work [27 \ \26 \ [28] has been the study of the ratios of partition functions 
(II. 6p and (II. 8p at large distances r = \xi — X2I (the infrared (IR) region) for 1-l-l-dimensional 
integrable QFTs. Integrability means that in these models there is no particle production in 
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any scattering process and that the scattering (S) matrix factorizes into products of two-particle 
S'-matrices which can be calculated exactly (for reviews see e.g. [30l EU [32l [33l [Mj ) . This is the 
factorised-scattering theory of integrable models. Since the scattering matrix and the particle 
spectrum fully encodes the local definition of QFT, it is also possible to incorporate the presence 
of boundaries in an integrable model defined in the factorised-scattering way. The study of 
integrable QFTs with boundaries has attracted a lot of attention in the last two decades (see 
e.g. [351 EH EZl [291 EH [39] ) . In our work on the boundary case we have made extensive use of the 
results of Ghoshal and Zamolodchikov particularly the explicit realization of the boundary 
state which they proposed. 

Taking the known S-matrix of a model as input it is possible to compute the matrix elements 
of local operators (also called form factors) . This is done by solving a set of consistency equations 
[40\ I41j . also known as the form factor bootstrap program for integrable QFTs. In [27J, this 
program was used and generalised in order to compute (11.60 in the case of integrable models with 
diagonal scattering matrix (that is, without backscattering) . In order to do this, the two-point 
function in (II. 6p was expressed as a sum in terms of form factors of the twist fields involved (an 
expansion using a decomposition in energy-momentum eigenstates). This was then extended 
to models with backscattering, such as the sine-Gordon model [26], to integrable models with 
boundaries in [28], and some aspects were generalised to non-integrable models [42]. 

One of the most interesting results of our works [27t [^ WI\ has been the identification of 
the next-to-leading order correction to the large-distance (large-r) behaviour of the entropy of 
all (unitary) massive two-dimensional theories, that is, the third term in the following large-r 
expansion: 

ST'^ir) = -| log(em) + U'^""'' - - ^ i^o(2rm„) + O (e'^^"^) . (1.12) 

Here, nia are the masses of the £ particles in the QFT model, with m = mi < rua, and [/™°dei 
the model-dependent constant introduced above. The first two terms are the expected saturation 
of the entanglement entropy, but the interesting feature is the universal third term, where we 
see that the leading exponential corrections are independent of the scattering matrix, and only 
depend on the particle spectrum of the model. This is quite striking: for instance, a model 
of i free particles of masses rUa will give the same leading exponential corrections as one with 
interacting particles of the same masses. The result ()1.12p was first obtained using integrable 
QFT methods [27^ I26j. then, even more strikingly, it was understood to hold as well outside of 
integrabilty [12]. 

Infrared corrections were also studied for integrable models with boundary j28j , in which case 
they are always model-dependent, in particular through the reflection matrices off the boundary. 

The results described so far have been extended further for the particular case of the Ising 
model. For this model the particular form of all infrared corrections to the entropy with and 
without boundary has been obtained in [28j. In fact, for this model without boundary, like for 
other free-field QFT models, there is an alternative powerful way of studying the entanglement 
entropy, see the review [33] . 

This paper is organized as follows. In section 2 we review the relationship between partition 
functions on multi-sheeted Riemann surfaces, correlation functions of branch-point twist fields 
and the entanglement entropy. Employing these relationships, we provide expressions for the 
bi-partite entanglement entropy of H-l-dimensional quantum field theories, both in the bulk and 
boundary cases. In section 3 we introduce the form factor program for branch-point twist fields. 
We explain how it can be employed to obtain the form factors for these twist fields in integrable 
models, and how these form factors can be checked for consistency against conformal field theory 
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results. For the two-particle form factors we generalize this program also to 1+1-dimensional 
non-integrable QFTs. We identify the general structure of the two particle form factors and, in 
integrable cases, suggest how higher particle form factors may be obtained from lower particle 
ones. For the Ising model, we give closed formulae for all higher particle form factors. In general, 
we recall how the form factors can be regarded as building blocks for correlation functions. The 
correlation functions can then be expressed as series where the leading contributions at large 
distances arise from the lower particle form factors. In section 4, we used these form factor series 
in order to analyse the entanglement entropy, both in the bulk and boundary cases. In the bulk 
case, our most important result is the universal expression ()1.12p for the next-to-leading order 
correction to the entanglement entropy at long distances. Both for the bulk and boundary cases 
of the Ising model, we evaluate all higher order infrared corrections to the entanglement entropy. 
In the boundary case, we find the precise relationship between the ultraviolet leading behaviour 
of the entanglement entropy and the boundary entropy introduced by Affleck and Ludwig |44j . 
Finally, in section 5 we summarise our conclusions and outlook. 

2 Replica trick and entanglement entropy 

2.1 Partition functions on multi-sheeted Riemann surfaces 

From the considerations in the introduction, it is clear that for the study of the entanglement 
entropy in the scaling limit, we must study partition functions of (euclidean-signature) quantum 
field theory on multi-sheeted Riemann surfaces. We will come back to the precise relation in 
the next section, but for now, let us discuss such partition functions. In particular, we wish to 
introduce the concept of branch-point twist fields, following |27j . 

The partition function of a model of 1+1-dimensional QFT, with local lagrangian density 
£[93](x), on a Riemann surface TZ is formally obtained by the path integral 



Here, [d(p]'ji is an infinite measure on the set of configurations of some field (p living on the 
Riemann surface TZ and on which the lagrangian density depends in a local way, and xti is a point 
on the Riemann surface. Consider Riemann surfaces with zero curvature everywhere except at a 
finite number of points; the points where the curvature is non-zero are branch points. In the case 
where the initial quantum system we consider is on the line, the Riemann surface forms a multiple 
covering of M^; in the case where the model is on the half-line, it forms a multiple covering of 
the half-plane, which we will take to be the right half-plane, = {(x, y) | x > 0, y G M}. For 
simplicity, we will only consider the case in the following discussion, but the case is 
entirely similar. Since the lagrangian density does not depend explicitly on the Riemann surface 
as a consequence of its locality (and the fact that the curvature is zero almost everywhere) , it is 
expected that this partition function can be expressed as an object calculated from a model on 
M?, where the structure of the Riemann surface is implemented through appropriate boundary 
conditions around the points with non-zero curvature. Consider for instance the simple Riemann 
surface M.n,xi,x2 with xj = (xj,0), composed of n sheets sequentially joined to each other on 
the segment x E [xi,X2], y = (see Fig. [2]). We would expect that the associated partition 
function involves certain "fields"Lj at the points on given by xi and X2- These fields would 

*Here, the term "field" is taken in its most general QFT sense: it is an object of which correlation functions - 
multi-linear maps - can be evaluated, and which depends on a position in space - parameters x, y that transform 
like coordinates under translation symmetries. 
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be implementing the appropriate non-zero curvature at the branch points. 

The expression (12. 1|) for the partition function essentially defines these fields: it gives their 
correlation functions, up to a normalisation independent of their positions; we only have to put 
as many branch points as we have such fields in the correlation function. But if we insist in 
understanding this as the initial model on M?, this definition makes these special fields non- 
local. Locality of a field (used here in its most fundamental sense) means that as an observable 
in the quantum theory, it is quantum mechanically independent of the energy density at space- 
like distances. In the associated euclidean field theory, this means that correlation functions 
involving these fields and the energy density are, as functions of the position of the energy 
density, defined on and continuous except at the positions of the fields. The energy density 
is simply obtained from the lagrangian density. But clearly, bringing the lagrangian density 
around a branch point changes the value of the correlation function, since it gets to a different 
Riemann sheet. Hence, the fields defined by ()2.ip seen as fields in a model on with lagrangian 
density £[(/9](x) are non-local: the lagrangian density is not well-defined on E?. Locality is at 
the basis of most of the results in QFT, so it is important to recover it. 

In order to correct the problem, note that if we defined a new lagrangian density on M? at 
the point (x, y) by simply summing the initial lagrangian density '^'j=i over all the 

points of the Riemann surface TZ = M.n,xi,x2 that project onto (x, y), then with respect to 
this new lagrangian density, the fields defined above would be local. Hence, the idea is simply to 
consider a larger model: a model formed by n independent copies of the original model, where 
n is the number of Riemann sheets necessary to describe the Riemann surface by coordinates 
on M^. Let us take again the simple example of M.n,xx,x2- We re-write (12. ip as 

Z„(xi, X2) ^ Mn,x,,x,] = f [d^i ■ ■ ■ dy,J^2 e" '^'^ (^Kl (2.2) 

JC(xi;x2) 

where C(xi,X2) are continuity conditions on the fields ip^, . . . ,ip^ restricting the path integral: 

C(xi;x2) : V9j(x,0+) = v'j+i(x,0") , x e [xi, X2], z = 1, . . . , n, (2.3) 

where we identify n + i = i. What appears in the action of that path integral can now be seen 
as the lagrangian density of the multi-copy model, 

, . . . , ip^] (x, y) = (x, y) + . . . + C[ip^] (x, y) . 

The energy density of that model is the sum of the energy densities of the n individual copies. 
The expression (|2.2p with (j2.3p does indeed define the insertion of local fields at xi and X2 
in the multi-copy model, since the energy density is the same on both sides of the segment 
X G [xi,X2], y = according to the conditions C(xi,X2). 

The local fields defined by (j2.3p are examples of "twist fields". Twist fields exist in a QFT 
model whenever there is a global internal symmetry a (a symmetry that acts the same way 
everywhere in space, and that does not change the positions of fields): dxdy C[aip]{x,y) = 
J^2 dxdy C[ip]{x,Y) for some lagrangian density C[ip]{x). Their correlation functions can be 
formally defined through the path integral: 



(2.4) 



(T^(x) • • •)^jg2 oc / [dip]^2exp - dxdy C[ip]{x,y) 

JCa{x) I 

where • • • represent insertions of other local fields at different positions. The path integral 
continuity conditions produce a cut on a half-line starting at the point x = (x, y): 

Ca{x) : (/?(x',y+) = cr(/?(x',y") , x' G [x, cx)) . (2.5) 
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The proportionality constant is an infinite constant that is independent of the position x and of 
those of the other local fields inserted, present in order to render the path integral finite (so that 
it may represent a correlation function). For insertion of many twist fields, we just add more 
continuity conditions on different half-lines, starting at different points. The fact that cj is a 
symmetry ensures that 7^ is local, since it ensures that the energy density is continuous through 
the cut produced by the continuity condition. Also, it ensures that the result is in fact invariant 
under continuous changes of the shape of this cut, up to symmetry transformations a of the 
local fields that are being swept. This is because inside a loop, we can always apply a symmetry 
transformation without changing the result of the path integral, up to transformations of local 
fields present inside the loop. In this way, we can modify the continuity conditions through the 
loop. By drawing a line, where the usual continuity holds, starting and ending on a twist-field 
cut (or possibly on its end-points), we form a loop with the twist-field cut. Applying a symmetry 
transformation inside this loop, we erase part (or all) of the twist-field cut and make the line a 
new twist-field cut. In this way, it is possible to move the twist-field cut to any shape. Hence 
our choice of a half- line extending to the right in the definition of the twist field is just for 
convenience. 

A consequence of the formal definition ()2.4p is that correlation functions (7^(x)0(y) • • •)£ ]g2 
with some local fields 0{y) are defined, as functions of y (continuous except at positions of other 
local fields), on a multi-sheeted covering of M? with a branch point at y = x, whenever aO ^ O. 
More precisely, twist fields have the property that a clockwise continuous displacement of a local 
field 0{y) around x back to its initial projected point on is equivalent to the replacement 
O ^ aO in any correlation function. If aO ^ O, then O is said to be "semi- local" with respect 
to 7^. This twist-field property is satisfied by a large family of fields, not only the one obtained 
through the formal definition (12. 4p . For instance, we could have inserted a field (/^(a;) in the path 
integral, leading to the same twist property (this is a descendent twist field). However, the twist 
property, along with the condition that 7^ has the lowest scaling dimension and be invariant 
under all symmetries of the model that commute with a (that is, that it be a primary field in 
the language of conformal field theory), uniquely fixes the field 7^ up to a normalisation. These 
conditions lead to a definition that is in agreement with the path-integral definition ()2.4p . In 
fact, these conditions constitute a more fundamental way of defining the primary twist field than 
the path integral, as they do not require the existence of a lagrangian density. In particular, 
they lead to unambiguous definitions in any quantisation scheme (we will discuss what the twist 
condition is in the quantisation on the line when we discuss form factors). We will take this 
general point of view in the following, but we will continue to think of a model of QFT through 
its lagrangian density C for clarity. 

Twist fields associated to internal symmetries have been largely studied in the context of 
CFT: they correspond to twisted modules for vertex operator algebras [46 ^ 147 ^ 118]. at the basis of 
so-called orbifold models. In the context of massive integrable QFT, only the simplest (standard) 
cases are well known. The Ising order and disorder fields are Z2 twist fields in the free massive 
Major ana fermion theory, and in the equivalence between the massive Thirring model and the 
sine-Gordon model, the bosonic exponential fields of the latter are U{1) twist fields of the former. 
But in fact, for fields with more general twist propertes, the form factor equations of integrable 
QFT were written in |49j . The first extensive study of form factors of non-standard twist fields 
(see below) was done in \27\ I26j . and the first time twist fields form factors were considered 
beyond integrability was in [l2]; we will describe these works in the next section. 

In the n-copy model with lagrangian £ = there is a symmetry under exchange of 

the copies. The twist fields defined by (|2.2p . which we call branch-point twist fields, are twist 
fields associated to the two opposite cyclic permutation symmetries i 1— > i + 1 and i + 1 t-^ i 
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Figure 3: The effect of T on other local fields. 



{i = 1, . . . ,n, n + 1 = 1). We will denote them simply by T and T, respectively: 

T = %j , a : i i + 1 mod n 

T = 7^-1 , : i + 1 i modn 

(see Fig. [3] for the case T). More precisely, we have 

Zn{xi,X2) OC {T{xi)f{x2))^(u)^^2. (2.6) 

This can be seen to be correct by observing that for x G [xi, X2], consecutive copies are connected 
through y = due to the presence of the cut produced by T(xi), whereas for x > X2, the 
additional cut produced by T(x2, 0) cancels this, and copies are connected to themselves through 
y = 0. Hence, indeed only a finite cut between xi and X2 remains. The precise proportionality 
constant was discussed around equation p.6|) . 

More generally, the identification holds for correlation functions in the model C on Ain,xi,x2i 
this time with an equality sign: 

/m 1. f •^ \ (T(xi)t(x2)0i(?/)---)^(n)^R2 

{0{yTi on sheet i) ■ ■ -jcMn..^,.^ = (2-7) 

{T{xi)T{X2)}c(n),R2 

where Oi is the field in the model C^""^ coming from the z*^ copy of C, and y is the projection of 
y-n onto M^. 

Note that it is easy to transform the twist field T into T, and vice versa: we only have to 
apply the "fiip" symmetry transformation by which the order of the copies is inverted (this is 
another element of the permutation symmetry group). Note also that this construction can also 
be generalised to Riemann surfaces with more branch points is straightforward, but this will not 
be needed here. 

The conformal dimension of branch-point twist fields is an important characteristic of these 
fields. It was essentially calculated in |22j . although branch-point twist fields were not in- 
troduced; only the non-local fields discussed above were considered. Here, we reproduce the 
derivation of [27], which makes explicit reference to the branch-point twist fields, but otherwise 
follows closely |22j . 

Consider the model to be a conformal field theory (CFT). Then also is a CFT. There 
are n fields Tj{w) in >C^") that correspond to the holomorphic stress-energy tensors of the n 
copies of C, and in particular the sum T^"'\w) = X]j=i^i(^) is the holomorphic stress-energy 
tensor of The central charge of C^"'^ is nc, if c is that of C. 

Consider the holomorphic stress-energy tensor T{w) in C We can evaluate the one-point 
function {T{w))c,Mn,xj^,x2 making a conformal transformation from z in to w in A4n,ai,a2 
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(here z and w are complex coordinates, and Wj = Xj + iyj) given by 

' w — W\ 



We have 



W — W2 



where the Schwarzian derivative is 

{z,w} 

Using {T{z))£^^2 = 0, we obtain 
{T{w))c,M 



(dz \ c 
^1 {T{z))c,-ig2 + —{z,w} 

z"V-(3/2)(z")2 



'^2 



c{'n? — 1) (-wi — 

2Aii? (w — wiY^w — 102)"^ 



Since, by (|2.7|) . this is equal to {T{xi)T{x2)Tj{w)) ^(n) ^2 / {T {xi)T {X2)) /^(v.) ^2 for all j, we can 
evaluate the correlation function involving the stress-energy tensor of C^'^'^^ by multiplying by n: 

c{n'^ — 1) {wi — W2)'^ 



{T{xi)f{x2))c(n)^R2 24n {w -wi)'^{w -102^ 

From the usual CFT formula for insertion of a stress-energy tensor 



{T{x^)f{x2)T(-Hw))^, 



If — (it; — wiY w — W2 dw2 {w — W2Y'' 

we identify the scaling dimension of the primary fields T and T (they have the same scaling 
dimension) using {T {xi)'T {X2)) /^in) ^^2 = \wi —W2\~'^'^'^: 

<i„ = ^(n-i). (2^8) 

This dimension is the lowest possible dimension for fields with the branch-point twist property, 
as it is the dimension of a primary field in this family. 

Hence, the branch-point twist fields T and T are unambiguously defined by specifying their 
associated symmetry, a and ct~^, namely that they are invariant under other global symmetry 
transformations that commute with cr, and by specifying their scaling dimension to be d„ (12. Sp . 
This is a definition up to normalisation; we will come back to the normalisation in the next 
section. There are of course many other twist fields associated to other elements of the permu- 
tation symmetry group, and these may have smaller scaling dimension than but these are 
not the branch-point twist fields that are used in order to represent the partition function on 
multi-sheeted Riemann surfaces. 
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2.2 Bulk entanglement entropy 

Partition functions on Riemann surfaces with branch points can be used in order to evaluate the 
entanglement entropy. We start by considering the bulk case, where the quantum model is on 
the line, and the corresponding (euclidean) QFT model is on M^. We are ultimately interested 
in the scaling limit of the ground-state entanglement entropy, (|1.2|) . However, we will provide 
here general arguments valid not only for the ground state, by also for excited states. Hence, 
in place of |gs) in (|1.2|) . we will take some arbitrary, possibly excited state \ip). In the scaling 
limit, the ground state |gs) maps to the QFT vacuum state |0), but excited states map to QFT 
asymptotic states characterised by the particle content and their momenta (in general, we will 
keep the notation 1-0) for the corresponding state in QFT). We will later specialise \tp) to the 
ground state. 

The main idea in order to evaluate the entanglement entropy in the scaling limit of a quantum 
model is to use the replica trick [20\ [2H [22} ; we will follow more precisely the method of |22j . 
That is, we use the identity (|1.4|) . where we start by considering the trace for positive integer n, 
evaluate it in QFT (in the scaling limit), then take the appropriate analytic continuation in n 
in order to evaluate the limit of the derivative. Considering positive integer n in QFT is useful, 
because there, the trace is directly related to a partition function on a multi-sheeted Riemann 
surface, studied in the previous subsection. The way this works is as follows. 

First, in order to construct in the scaling limit, consider the QFT Hilbert space as a space 
of field configurations {ip{x), x E R} on R. A state can be written as a linear combination of 
field configurations on R where the coefficients are obtained by path integrals on the lower half 
of the R^ plane, R^ = {(x, y) | x G R, y < 0}. The boundary condition on y = is determined 
by the coefficient we are looking at, and the asymptotic condition y — > is determined by I?/'). 
That is, 

[#]m \mm, m) = ^ Jmui e''"""^^^ (2.9) 

where ^^[v?] = (fx C[ip]{x) is the action of the model on Ai, and the condition C is 

^ . r (^(x,0)=(/.(x),xGM 

W(x,y ^ -00) ~ /^(x,y), X G ] 

Here, /^(x, y) represents the asymptotic condition corresponding to the state \ip); for instance, 
/^(x, y) = for the vacuum state, otherwise it reproduces wave packets corresponding to asymp- 
totic particles. The number Zi is the "partition function" of the QFT model in the state \^): 



c 



where 




-00) ~ /v,(x,y), X G 
oo)~/;(x,y),xGl 



It is a true partition function if lip) is the ground state; otherwise, it corresponds to an excited- 
state normalisation. The factor \l\fZ\ is inserted here in order that we maintain the proper 
normalisation (V'lV') = 1 (the states |0) are taken with the natural delta- functional normalisa- 
tion). Note that the coefficients in the expansion of the dual vector are simply obtained 
from (12. 9p by complex conjugation, along with the inversion of the imaginary time y ^ — y - 
this naturally gives a path integral on the other half of . 
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Next, the density matrix is constructed by tracing \tp){ilj\ over degrees of freedom on 
A. We sum the diagonal matrix elements of but only in the space A, formed by field 

configurations on A. This has the effect of "connecting", on the part A of their boundaries, the 
two surfaces on which the path integrals used to represent matrix elements of {ip) and of are 
defined. Hence, we get a path integral on the whole M^, with continuity on A, but with an open 
slit on A. We are left with a matrix element of p^, determined by field configurations on both 
sides of the slit on A. That is, we find 

{^aIPaM = J [#a]a {cPAAA\i^)mA,'P'A) = j^jt^¥']R2\(A,0) e'^^^VC^.o) (2.10) 
where 

C" 



99(x,0-) = 0^(x), x£A; ^(x,0+) = 0^(x), x G A 
V{^,y ^ -oo) ~ /v,(x,y), X e M; 99(x,y ^ oo) ~ /^(x,y), x E 



The n^^ power of 



Pa, 



{<Pa\Pa\<P'a) = J [#1#2 • • • d(l)n-l]A{(pA\pA\(Pl){<Pl\PA\<p2) ' ' ' IPaI0a) > 

can be obtained by taking n copies of with their own independent path integrals over ipj, j = 
1, . . . , n, and by connecting in a sequential way one side of the slit on A of one copy to the other 
side of the slit on A of the next copy (that is, having continuity of the fields ipj through the slit 
A in this sequential way). Finally, taking the trace connects the last copy to the first, so that 
we obtain the partition function on a multi-sheeted Riemann surface. Taking A to consist of 
only one interval for simplicity, with end-points at positions xi and X2 in M?, we then find, in 
agreement with (II. Sp . 

TrAPl = (2.11) 

where the partition function is (j2.2p . In this derivation, we referred to the definition of the QFT 
model given by its lagrangian density C, but as mentioned in the previous subsection, this is 
just for clarity of exposition; the existence or not of a lagrangian density does not affect any of 
the results. 

As we saw in the previous subsection, the partition function above can be computed as a 
two-point correlation function of local fields in the model given by the lagrangian density 
using ()2.6p . In the case of the ground state, using (II. 6p . we have 

ST^'ir) = - hm -^Z„e2rfn(o|T(xi)t(x2)|0). (2.12) 



In the case where \ip) is an excited state, the relation (j2.6p still holds, because the derivation of 
the previous subsection did not make any reference to the asymptotic conditions on the field ip as 
y iboo. However, in (j2.6p . we have to understand the correlation function as an excited-state 
diagonal matrix element of the product of twist fields, in order to implement the appropriate 
asymptotic conditions in the path integral. Hence, what replaces (ll.6p is 



so that we obtain 



^H^^^ = Z^e''-{mx^)f{x2m. (2.13) 



ST\r) = - lim ^Zne^^-{^p\r{xl)f{x2m. (2.14) 
n^i an 
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In this formula, the normahsation of is not the standard one in the case of an excited 
state, because of the normahsation factor in (I2.13|) . Indeed, this factor normahses away 

the infinite-volume divergencies coming from colliding rapidities of asymptotic states in the 
disconnected terms of the two-point function. We will not discuss further here the case of 
excited states, and concentrate solely on the ground state entanglement entropy. 



2.3 From bulk to boundary entanglement entropy 

We are now interested in the situation where the system is on the half-line x > and the 
connected region A has an end-point on the boundary of the system x = 0. However, it will 
be most easy and instructive to start with a slightly different problem, where the region A lies 
entirely in the bulk, with two boundary points on the half line x > 0. Since the system is on 
the half-line, one needs to provide an additional boundary condition at x = in order to fully 
define the model. There are various ways of implementing such a boundary condition. From 
the viewpoint of the path integral, the boundary condition is implemented by a restriction on 
the allowed values of the fields <^(x) and its derivatives at the boundary x = 0, along with, 
possibly, an extra term Sb in the action that is supported on x = 0, 5*6 = / dyCB['^]{j)- From 
the viewpoint of the quantised theory on the half-line x > 0, with y G M the imaginary time, 
the boundary condition determines the whole Hilbert space, the vacuum |0)s and all excited 
states. Finally, crossing symmetry gives, from the latter, the viewpoint of the quantised theory 
on the full line y G M with imaginary time x > 0. There, the boundary condition corresponds 
to a boundary state |-B), a state in the usual Hilbert space on the full line. 

The derivations of the previous two subsections, connecting the entanglement entropy to 
partition functions on multi-sheeted Riemann surfaces, and then connecting the latter to correla- 
tion functions of branch-point twist fields, can be directly generalised to the boundary situation. 
First, retracing the steps of the last subsection in the path- integral formulation of the model on 
the boundary, we obtain again (I2.11|) . where now the partition functions are path integrals over 
configurations of the field (p{x) with x G M^, with a boundary condition at x = and possibly 
an extra boundary term Sb in the action. Second, the derivation of subsection 12.11 can also be 
directly reproduced, and we obtain, instead of (jl.6p . 



= Z„e2^"B(0|r(xi)r(x2)|0)B. (2.15) 

Here |0)b is the ground state in the n-copy QFT model on the half line, and T and T are 
the twist fields with the same fundamental definition as before (the twist property, invariance 
under other symmetry transformations, and lowest scaling dimension), but as operators on the 
half-line Hilbert space. Hence we obtain 

^boundary (^^^^^^ = - lim ^ Z^E^'^- B^ {^l)^ {X2m B ■ (2.16) 

n-»i an 

From this, there are two ways to obtain the entanglement entropy S^^^'^^'^^ {r) for a region 
A starting from x = and ending at a distance r from it. First, we could consider the limit 
xi — > (0,0) in S^^^^'^^'^^ {xi,x) with x = (r, 0), making the region A approach the boundary. 
As the twist field at xi = (xi,0) approaches the boundary, the correlation function diverges, 
because the presence of the boundary changes the regularisation necessary around the branch 
point. A way to evaluate the divergency is to use boundary conformal field theory, which applies 
in massive models when a local field is near to a boundary. It tells us that for small xi there 
is a power law determined by the conformal dimension of T: b(0| •••T(xi)|0)b oc x^*^" |50j . 
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We may then define T(0)\0)b as limx^-^o xf"T(xi)|0)B- This appropriately regularised operator 
T(0) is simply proportional to the unitary operator performing a cyclic Z„ transformation, 
since its branch cut, through which Z„ transformations are performed, now extends through 
the whole space. But since \0)b is invariant under such a transformation, the action of T(0) on 
\0)b gives \0)b back. Moreover, since in the definition of T(0) we already took into account a 
renormalisation of the field, there is no regularisation e'^" associated to that field in the relation 
between the partition function and the two-point function. Hence, we find, with appropriate 
choice of proportionality constants, 

^boundary ^ _ J- Z^s''- b{0\T {x)\0) B ■ (2.17) 

n^i dn 

The power of e guarantees that the scaling dimension of the quantity that is differentiated is 0. 

It is important to note that the result of the limit xi (0, 0) is a function only of the 
distance r between x and the boundary. Indeed, T(0) is a unitary operator preserving \0)b, 
hence can be put at any imaginary time without changing the result. 

Second, we may take the limit X2 — > oo in S^^^{x,X2), with again x = (r, 0). This should 
lead to the same result for the entanglement entropy, because of the symmetry Sa = S^, up 
to an additive term corresponding to the contribution to the entanglement entropy around the 
boundary point at oo. In the limit X2 — > oo, the two-point function in (|2.12|) reduces to its 
disconnected part: 

B(0|T(x)t(x2)|0)B ~ B(0|T(x)|0)BB(0|t(oo)|0)B . (2.18) 

In the second factor, the twist field does not feel the presence of the boundary, hence this 
expectation value can be replaced by its expectation value in the model without boundary, 
(0|T|0). Dividing out this factor (which corresponds to subtracting the contribution to the 
entanglement entropy around the point at oo), fiipping the sheets in order to transform T into 
T, and using the appropriate branch-point regularisation, we find again (I2.17p . 

Some of these considerations, and in particular the calculations that we will present below, 
are made clearer by using crossing symmetry: taking the half-line x > to be imaginary time, 
and the full line y G M the space direction. In this picture, as we mentioned, the boundary 
condition is implemented as a boundary state on the Hilbert space on the full line: 

B(0|T(xi)t(x2)|0)B = (0|T(xi)t(x2)|i?). (2.19) 

The boundary state |-B) is in the past at imaginary time x = in the Hilbert space of the model 
on the full line, and the twist fields are placed at imaginary times xi and X2. More precisely, 
the boundary state is the n-fold tensor product of single-copy boundary states. The state (0| is 
again the ground state of the n-copy model on the line, corresponding to asymptotic vanishing 
fields conditions at positive infinite times. No factor occurs in using crossing symmetry since the 
branch-point twist fields are spinless. The normalisation of the boundary state \B) is such that 
{0\B) = 1. By any of the two ways explained above, we obtain for the entanglement entropy: 

^boundary = - lim Z^s'^" {0\T {x)\B) . (2.20) 

where again the distance between the point x and the boundary is r. 
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2.4 Normalisations of branch-point twist fields and boundary states 



In order to obtain a universal result for the entanglement entropy, we need to fix the short- 
distance behaviour as in (|1.9|) in the bulk case, and fix the large-distance behaviour as in 
in the boundary case. As was explained in the introduction, this provides a QFT definition for 
£ (in general, different for the bulk and boundary cases), hence a universal definition for the 
finite terms in Sa- 

The general definition of the branch-point twist fields as explained in subsection 12.11 did not 
say anything about the normalisation of the fields. This normalisation does not depend on the 
positions of the fields, but may, and in general does, depend on the number of copies n - that is, 
a change of normalisation may change the n-dependent number Z„ in ()2.12p and ()2.20p . Let us 
choose the CFT normalisation: this is the normalisation by which the short-distance behaviour 
of the two-point function is exactly given by the CFT two-point function, with coefficient 1: 

(0|T(xi)T(a;2)|0) ~ |xi -xar^"^" as m|xi - X2I ^ 0. (2.21) 

This two-point function just depends on the euclidean distance |xi — 3:2 1, thanks to euclidean 
translation and rotation invariance (or Poincare invariance in Minkowsky space) , and to the fact 
that the twist fields are spinless. With this normalisation, and with a fixed relation between e 
and the correlation length ^, we have a fixed function Z„ in (jl.6p . In particular, Zi = 1, because 
with the CFT normalisation, the twist fields just become the identity operator 1 at n = 1. We 
may change the definition of e by rescaling it in an n-independent way, so that we change its 
relation to the correlation length, but keep it independent from n. This rescaling is absorbed 
into Zn, and we may do so in order that dZn/dn = at n = 1, as said in the introduction. With 
this rescaling of e, it is easy to see that taking the derivative in (j2.12p with the short-distance 
behaviour (I2.21|) . we immediately obtain the UV asymptotic given by (II. 9p for the entanglement 
entropy. 

Concerning the boundary case, we need to choose the normalisation of the boundary state \B) 
in order that with the CFT normalisation (j2.2ip . the expression (j2.20p gives rise to the correct 
IR asymptotic in (II. lip . In that IR asymptotic, there is the constant that appears, a 

constant that is defined in the bulk situation. At large r, the two-point function in the bulk 
factorises into a product of one-point functions (this is asymptotic factorisation of correlation 
functions of local fields): 

(0|T(xi)T'(x2)|0) ~ ((0|T|0))2 as m|xi-X2|^oo (2.22) 

where we used the fact that T is obtained from T by a permutation symmetry transformation 
that keeps the vacuum state invariant. Using this into ()2.12p . we find the correct IR behaviour 
of ([H]), with 

^modei ^ _2 lij^ A. ( m-'^n /olTIO)) . (2.23) 
n^i an \ J 

Hence, this constant is related to the one-point function of twist fields in the CFT normalisation. 
Returning to the boundary case, with the choice of boundary state normalisation 

(0|B) = 1, (2.24) 

we see that a decomposition into energy and momentum eigenstates, keeping the ground state 
only, gives 

(0|r(x)|5) ~ (0|r|0) as mr\ 00 (2.25) 
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where r is the distance between x and the boundary. Hence, we indeed find, from (j2.20p . the 
correct IR behaviour given by (II. lip for the entanglement entropy in the boundary case. 

With this normalisation, the model-dependent and boundary-condition-dependent constant 
V{k) is then universal. Since it occurs at small distances, it is derived from the UV behaviour of 
the one-point function in (j2.20p . This behaviour is itself completely controlled by the conformal 
UV fixed point, hence at short distances that one-point function is a CFT one-point function. 
This means that the constant V{k) only depends on the conformal boundary condition to which 
the system flows as mr — > 0. This was seen explicitly in the Ising model in [5H], which will be 
reviewed in section |H 

3 The form factor program for branch-point twist fields 

We will now describe how to evaluate correlation functions of branch-point twist fields, involved 
in the universal part of the entanglement entropy in the bulk and boundary cases through the 
expressions (|2.12p and (|2.2Up . The method we will use is that of form factor expansion; this is 
a representation of the two-point function (in the bulk case) or the one-point function (in the 
boundary case) that uses the decomposition of the identity operator into projection operators 
on energy and momentum eigenstates. It is an effective large-distance (large-r) expansion, and 
in principle, its re-summation gives a representation valid at all non-zero distances. 

3.1 Form factor expansion in l-|-l-dimensional massive QFT 

We turn to the description of QFT on Minkowski space-time in terms of its Hilbert space of 
asymptotic relativistic particles. This will allow us to introduce the method we use in order to 
evaluate the correlation functions of branch-point twist fields: the form factor expansion. 

In the context of 1 -|- 1-dimensional QFT, form factors are defined as tensor valued functions 
representing matrix elements of some local operator 0{x) located at the origin x = between a 
multi-particle in-state and the vacuum: 

Ff l^-^^(^i, ...,6k) := {0\Ome^, . . • , ^fc>J^,..„^^ ■ (3.1) 

Here l^i,...,^^)™ represent the physical "in" asymptotic states of massive QFT. They 
carry indices /Xj, which are quantum numbers characterizing the various particle species, and 
depend on the real parameters 6i, which are the associated rapidities. The form factors are 
defined for all rapidities by analytically continuing from some ordering of the rapidities; a fixed 
ordering provides a complete basis of states, for instance 9i > 62 > ■ ■ ■ > 0^. Form factors 
also depend on the number of copies n, but for simplicity we will not write this dependence 
explicitly. Moreover, because of relativistic invariance and spinlessness of the branch-point twist 
fields, they in fact only depend on the rapidity differences; in the two-particle case, they then 
become functions of only one variable, 6 = 9i — 62. 

The twist property associated to the permutation symmetry in the general definition of 
branch-point twist fields was explained in subsection 12. 11 It is at the basis of the main properties 
of the form factors of branch-point twist fields, and ultimately of our large-distance result about 
the entanglement entropy. It implies the following fundamental exchange relations for branch- 
point twist fields as operators on the space of asymptotic states. If ^'i, . . . , are the fields 
associated to the fundamental particles of each copy of the original model, then the equal time 
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= y") exchange relations between T and ^'i, . . . , can be written in the following forrr0 

^i{y)T{x) = T{x)^^+i{y) x^ > y\ 

^iiy)T{x) = T{x)^,{y) x' < y\ (3.2) 

for i = 1, . . . ,n and where we identify the indices n + i = i. The twist field T has similar 
properties as T, with the difference that its exchange relations with the fundamental fields of 
the theory are given by 

'^i{y)f{x) = t(x)^',_i(y) x^>y\ 

^^{y)f{x) = f{x)^,{y) x^<y\ (3.3) 

instead of (j3.2|) . The connection with the previous section is obtained by recalling that in going 
from the Hilbert space description to the path integral description, the order of operators is 
translated into time-ordering (or y-ordering in euclidean space), and that left-most operators 
are later in time. 

For our present purpose, what is important is that these exchange relations, along with the 
uniqueness of the branch-point twist fields explained in subsection 12. H give us the following 
relation between branch-point twist fields, as operators on the space of asymptotic states: 

t = (3.4) 

In order to obtain large-distance expansions, we simply introduce a complete sum over quan- 
tum states 



/ '^^\2nf'' 1^1' • • • ' ^'^^Mi.-./^fe Mfe,...,/.i(6'fc, ■■■,di\ (3.5) 

'==l^i---^*6»i>6»2>->6»fc 

The two-point function in ()2.12p is then expressed as 



k 



(0|T(xOt(x2)|0)=5] 5] / ^'^;r ^ ' K0|T(0)|^i,...,^,),_,j2 

(3.6) 

where r = \xi — X2\ is a space-like distance between the points xi and X2 in Minkowsky space- 
time. Here, we used the fact that T is the hermitian conjugate of T ()3.4p . and we used translation 
covariance in space and time in order to extract the exponential factor containing all the position 
dependence. We also used relativistic invariance and spinlessness to say that the form factors 
only depend on the rapidity differences. In this way, we could do a simple shift of rapidities in 
the imaginary direction, so that only cosh^j's are left in the exponential, with coefficients the 
negative of the space-like relativistic distance times the mass of particle jjij. The integrals, then, 
become explicitly exponentially convergent (the large-rapidity asymptotics of the form factors 
is at most exponential in 0j's). In the same spirit, we could have simply thought of xi and X2 
as coordinates in 'M? (with imaginary time), used the euclidean rotation invariance in order to 
bring xi — X2 in the pure imaginary time direction (with a distance r), then used imaginary-time 
translation covariance to extract the real-exponential factor depending on r. The factor 1/(271)*^ 
in the expansion just tells us about the normalisation we have chosen for the asymptotic states. 



^Here we employ the standard notation in Minkowski space-time: with = 0, 1, with being the time 
coordinate and being the position coordinate. 
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Similarly, the one-point function in (j2.20p is expressed as 



{o\nx)\B) = E E 



d9i...d6k E n^Mj coshe^ 



X (O|T(O)|0i, . . . , 9k)t.,,...,,, A(^k, Oi\B)] . 



(3.7) 



Here, the easiest way to obtain this expansion is by considering that the distance r between 
X and the boundary is in the imaginary time direction, and using imaginary-time translation 
covariance in order to bring out the real-exponential factor. 

In the expansions ()3.6p and (13. 7p the terms with a larger number of particles are smaller 
because the lowest possible value of the argument of the exponential is —rkmi, k times the 
lowest mass mi. Moreover, as r increases, the larger-particle terms decrease faster; hence, these 
indeed are large-distance expansions. 



3.2 Form factor equations in integrable QFT 

We have found that (|2.2ip and (|2.24p provides the correct field and boundary state normali- 
sations in order to reproduce ()1.9p and (II. lip . We have also explained how the form factors 
can be used in order to calculate the correlation functions involved in the universal part of the 
entanglement entropy. Let us now provide more explanation as to the properties of form factors 
in the context of integrable models of QFT. 

The main characteristics of massive integrable models of QFT are that the number of par- 
ticles and their momenta set are conserved under scattering, and that the scattering matrix 
factorises into products of two-particle scattering matrices. The two-particle scattering matrix 
(or S'-matrix) is the solution of a set of consistency equations and analytic properties. These 
consistency equations and analytic properties are often strong enough, when combined with 
general properties of the model (symmetries, etc.), to completely fix the S-matrix in integrable 
models. In a similar way, the form factors can be completely fixed by solving a set of equations 
and analytic properties, which now depend on the two-particle S'-matrix. 

In this section we want to show how the standard form factor equations for 1 -|- 1-dimensional 
IQFTs must be modified for the branch-point twist fields. Here we will consider an integrable 
model consisting of n copies of a known integrable theory possessing a single particle spectrum 
and no bound states (such as the Ising and sinh-Gordon models). The interested reader may 
refer to 02] for the general case. We have therefore n particles, which we will denote by indices 
1, . . . ,n. The S'-matrix between particles i and j with rapidities 9i and 9j will be denoted by 
Sij{9i — 9j) (that it depends on the rapidity difference is a consequence of relativistic invariance). 
Particles of different copies do not interact with each other, so that the S-matrix of the model 
will be of the form 

Sij{9) = Si9)''^ V i,j = l,...,n, (3.8) 

where S{9) is the S-matrix of the single-copy integrable QFT. 

It is well known that exchange relations like ()3.2p and ()3.3p play an important role in the 
derivation of the consistency equations for the form factors. Generalising the standard arguments 
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Figure 4: A pictorial representation of the effect of adding 27ri to rapidity 9i in form factors of 
the twist field T. 

to the exchange relation ()3.2p . the form factor axioms are 

' ^'^'^'^ {. . . ,6i,9i+i, . . .) = Sf^^^.^_^{9ii+i)Fj^^ (. . . ,6'j+i,6'j, . . .j(3-9) 

^rK^2.../^fe(^^ + 27ri,...,efc) = Fjl^2-^"^i(02,...,^fc,0i), (3.10) 

Res F^l^^^'-'''' {00 + in, 00,6,..., Ok) = i F^''^-''^ {Ou ■ ■ ■ ,0k), (3.11) 



70 = f 



A: 

Res F^lf''^-''''{eo + i7r,eo,ei...,9k) = -illSi,^^{eo^)F^^'''■■■'''' {Oi, . . . ,9^). (3.12) 



70 = f 



i=l 



Here 9ij = 9i — 9j. Besides the simple poles whose residues are given in ()3.1ip and (j3.12p . 
the form factors are, as function of 9ij, analytic in the strip lm{9ij) G [0,27r) (there would 
be additional poles, with prescribed residues, if there were bound states; but we only consider 
integrable models without bound states here). The first axiom is in fact the same as for ordinary 
local fields. In the second equation, the crossing or locality relation, we introduced the symbols 
/ij = /ij + 1. As compared to the usual form factor equations, it is altered by the nature 
of the exchange relation and it now relates form factors associated to different particle sets. 
Finally, the last two equations generalise the standard kinematic residue equation to branch- 
point twist fields. Once more, the exchange relations ()3.2p are responsible for the splitting 
into two equations. Here, for later convenience, we wrote the equations in their general form 
valid also for many-particle models, where represents the anti-particle associated to fi. In the 
present case, the integrable model we started with has just one particle (so that fj, labels the 
copies) and therefore each particle is its own anti-particle. Pictorial explanations of the second 
and of the last two equations are given, respectively, in Figs. H] and \5\ 



3.3 Form factor equations in non-integrable 1+1-dimensional QFT 

Out of integrability, we do not have the luxury of the large simplification of the scattering 
amplitudes. However, two-particle form factors still satisfy certain properties. For simplicity, 
we will keep the notation as if there were only one particle type in the original model, so that 
fi represents the copy number, but attaching an extra index to does not affect any of the 
considerations or equations below. 
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Figure 5: The kinematic poles come from the structure of the wave function far from the local 
fields, at positive and negative infinity. Adding in to rapidity 6i puts the particle in the "out" 
region. With a particle in that region, there are delta-functions representing particles in the 
"in" region going through without interacting with the local fields. Those occur from the e*^^ 
form of the wave function at positive and negative infinity. But if the coefficients at both limits 
are different, S-e^^^ and S'+e*^^ with S- ^ S-^-, then there are also poles in addition to these 
delta-functions. Only these poles are seen in the analytic continuation i— > -|- itt. Different 
coefficients come from the semi-locality of the twist field and the non-free scattering matrix, as 
represented here. 

In the case of an ordinary, spinless local field O, the two-particle form factor ^^'^^^^^(01, 62) 
is a function of 9 = 9i — 62 with certain properties when considered as an analytic function of 
9, analytically continued from 9 > where it is the form factor with "in" asymptotic states as 
defined above. The properties of this function are: 

• the function analytically continued to negative rapidities 9 < along the real line is a form 
factor with an "out" asymptotic state: 

F^\^^'^^{e2 - 9,) mmOi,0^)Z,, ■ (3-13) 

• the function analytically continued by a 27ri shift from the real line gives: 

F^\^^^'^{9i -92 + 2ni) = Ff\^^^\92 - 9i) (3.14) 

• the function is analytic in the physical strip lm(0) G [0, 2ti) except possibly for poles on 
Re(0) = 0, lm{9) G (0, tt) U (tt, 27r) corresponding to bound states. 

Note that the physical strip in the rapidity plane, in the last point, is a double covering of the 
physical sheet in the Mandelstam's s-plane, obtained through s = rn\+m\+2mim2 cosh(0i— 6*2). 
This double covering has a symmetry under fiip with respect to the point ztt, as expressed by 
the second point above. For instance, for lm(0) = 2ti and Re(0) > 0, the function is describing 
a form factor with an "out" asymptotic state. 

The analyticity conditions above are somewhat strong conditions, but they are very far from 
allowing us to fix the two-particle form factors, because in general, the analytic structure outside 
of the physical strip is very complicated. However, in the case of the branch-point twist fields, 
the generalisation of these conditions are strong enough to give a result for the entanglement 
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entropy. The generalisation is obtained through arguments similar to those presented in the case 
of integrable models, whereby wave packets are exchanged or rotated around the point where the 
twist field lies, keeping only in mind that the exchange of wave packets does not lead to a simple 
two-particle S'-matrix in the non-integrable case. The first point above is unchanged by the 
twist property, as was equation (|3.9p in the case of integrable models, because it involves wave 
packet exchanges only on one copy of the replica model. Since in the replica model, particles on 
different copies do not interact even in general, non-integrable QFT, we still have here ()3.9p in 
the two-particle case when the copy numbers are different: 

pr\,.,^^_0^^pr\,^,2^0^ (^,^^2)- (3.15) 

The second point above is affected by the twist property, in exactly the same way as equation 
(|3.10|) was in the integrable case. Hence, equation (|3.1U|) still holds in the non-integrable, two- 
particle case: 

^r|/.iM2(^ + 27ri) = Fp^^^'{-6) (3.16) 

where we recall that fi = /u + 1. Finally, the third point above holds, up to one additional pole, 
the kinematic pole, as in the integrable case, equations (13. lip and ()3.12p . Indeed, the argument 
for this kinematic pole to be present, given in Fig. \5\ holds true when there is only two particles 
in the non-integrable case, because there is no exchange of wave packets necessary. Hence, we 
have 

Res Fp''''{e + i7r) = i(0|r|0), (3.17) 

e = o 

Res Fp''^{9 + i7r) = -i(0|T|0). (3.18) 
6 = 

Equations (IXTHIl . (IXTHIl . (IXTTD and (f3T8]) are what is left of ([Hj]), (l3ln]) . ([311]) and (l3l2D 
when integrability is taken away. As we will see below, in the integrable case, and only in this 
case, we can write exact form factors, but in general, the equations left in the non-integrable 
case are enough to obtain the exact result displayed in the introduction. 

3.4 Two-particle form factors 

Let us now analyse the equations for two-particle form factors in more detail. First, it is clear 
that by a cyclic exchange of the sheets, the branch-point twist fields are unchanged, so that form 
factors are invariant under such a cyclic exchange. Hence, it is sufficient to look at F^^''^{6), 
for instance. Second, we may repeatedly use (IXT^ and (nD3]l in order to obtain Fp^'^{0) for 
j = 2, 3, . . . , n. For instance: 

FP'\e + 27ri) = F^\'\-e) = F^'^ie), FP'\e + 2^i) = FP'\-e) = fP'\9), 

etc. Hence, we find 

FP'^'iO) = F^^'iff) V (3.19) 
F2^'\G) = Fl}^\e + 2^i{j-l)) i = 2,3,...,n. (3.20) 

Finally, the residue equations ()3.17p and (I3.18P as well as the analyticity conditions on the 
physical strip tell us about the analytic structure of -Fg {0) in the larger region lm(0) G 
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Figure 6: The structure of the function F2 (9) in the extended sheet Im(^) G [0,27rn], in 
the case n = 3. Crosses indicate the positions of the kinematic singularities. Shaded regions 
represent the physical sheets of the form factors F^^^^{9), F2^^'^{0) and F2^^^{0). 

[0,27rm), which we call the extended physical strip [27]. The function F2 (0) is analytic 
everywhere in that strip, except for simple poles on Re(0) =0: at 6 = iir, with residue i(0|T|0), 
and 6 = 2iTm — iir, with residue — i(0|T|0), and, possibly, on < lm{9) < vr and 27rn — vr < 
Im(6') < 27rn if there are bound states. This structure is depicted in Fig. [6] for the case n = 3 
and without bound states. 

This structure is all we can say about the two-particle form factors in the non-integrable 
case. In the integrable case, however, we can say much more; again, we concentrate here on 
diagonal scattering with a one-particle spectrum and without bound states for simplicity. As 
usual in the integrable context, we define the minimal form factors F^^^{6,n) to be solutions 
of the equations ()3.9p and (j3.10p for k = 2 without poles in the physical strip lm{6) £ [0,2it]. 
That is, 

FTiG) = FTi-mm = Flt^\^^^ -G) V J, k (3.21) 
where the S'-matrix is given by ()3.8p . Minimal form factors satisfy the same constraints (|3.19p 

_______ '7~| 11 

and (|3.20p . a consequence of which is that the minimal form factor F^^^ (0) must have no poles 
in the strip lm(9) £ [0, 27m] . From the equations above it is easy to deduce 

F^iO) = F2\-mO) = F2\-e + 2vrm). (3.22) 

These are the main equations from which we can determine all minimal form factors in integrable 
models. 

In order to develop a systematic procedure to solve these equations, it is useful to recall that, 
for a standard local operator the minimal form factor equations take the form 

fuiO) = hi{-d)S{ne) = M{-9 + 2TTi), (3.23) 

provided that the S'-matrix of the theory is given by S{n9). Thus given a solution to the previous 
equation, the function F^l^ {9,^) = fii{9/n) is automatically a solution of ()3.22p . 

In the context of integrable models, a systematic way of solving such type of equations has 
been developed whereby, given an integral representation for S{6), an integral representation 
of fii{0) can be readily obtained [30]. For diagonal theories, the integral representation of the 
S'-matrix takes the form 

S{9) = exp / -g{t)sm\i[ — 



(3.24) 
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where g{9) is a function which depends of the theory under consideration. A trivial consequence 
of the previous equation is 



S{n9) = exp 
and from here, it is easy to show that 
hi{9) =AAexp 



dt , , , . ^ fte\ 

I — qit/n) smh — 
Jo t ' ' WJ 



dt 



t sinh(nt) 



g{t) sin^ 



itn 



1 + 



vr 



where is a normahzation constant. Therefore, the desired solution is 



^Ii'n'W = /ii(W=AAexp 



dt 



tsinh(nt)' 



it 



g{t) sin 17 + — 



vr 



(3.25) 



(3.26) 



(3.27) 



In order to evaluate the full two-particle form factors, we must solve (|3.2ip with the poles in the 
extended physical strip mentioned before. The "minimal" solution, with in particular the least 
diverging behaviour at large 9, is 



(0) 



(T) sin 



,T\jk 



(0) 



2n sinh 



t7r(2(j-fc)-l)+6> 
2n 



sinh(MH(^^)^)F:if(zvr) 



(3.28) 



Note that the constant factor sinh(7r/n) guarantees that all form factors vanish for n = 1 as 
expected, since in that case the field T can be identified with the identity. 

For the field T, the exchange relations imply that form factors of the field T are equal to 
those of the field T up to the transformation i ^ n — i for each particle i. This means that 



T\(n~i){n~j) 



(9). 



This property can be combined with (|3.19|) - (|3.2U|) to show that 

-^T\ll, 



F^'^\9) 



\9). 



F;''\9) 



F^\^\9 + 2^i{j-l)). 



(3.29) 

(3.30) 
(3.31) 



3.4.1 Examples: the Ising and sinh-Gordon models 

The Ising and sinh-Gordon models are examples of the kind of integrable models for which we 
have specialised the results of this section: they have a single particle spectrum and no bound 
states. Their two-particle 5-matrices are. 



-1, 



and Ss.g{9) 



tanh\{9-i^) 



(3.32) 



respectively. In the sinh-Gordon model, the parameter B £ [0, 2] is the effective coupling 
constant which is related to the coupling constant (3 in the sinh-Gordon Lagrangian [51[ Wl\ as 



2/32 



87r + /3^ 



(3.33) 



under CFT normalization [53]. The 5-matrix [54^ [55l [56] is obviously invariant under the 
transformation B 2 — B, a symmetry which is also referred to as week-strong coupling 
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duality, as it corresponds to B{P) — > i?(87r/3~^) in (j3.33p . The point B = 1 is known as the 
self-dual point. 

Form factors of local fields of the Ising model were first computed in [10} [57] and later on in 
[58] for so-called descendant fields. A solution of ()3.2ip for j = A: = 1 is given by 

-^min^(^) ~ — ^sinh 5 fo'^ Ising field theory. (3.34) 

In accordance with our previous arguments, this is the standard minimal form factor already 
employed in 001 [57], with 9 6/n. 

Form factors of the sinh-Gordon model were first computed in [59j. The program was 
thereafter extended to other operators in [60] and more recently in [61j. The S'-matrix above 
admits an integral representation which is given by (I3.24p . with 



_ 8 sinh (f) sinh (i (l - |)) sinh (1) 



sinhi 

Therefore, the minimal form factor is given by 



(3.35) 



^Ii'n'W = exp 



°° dt sinh i§ sinh ^^^^ cosh t (n + 



for the sinh-Gordon model, 

Iq t sinh(nt) cosh g 

(3.36) 

where we have chosen the normalization J\f = F^^^ (zvrn). This integral representation can 
be turned into an infinite product of ratios of Gamma-functions as described in [27]. As a 
consistency check, it is quite easy to show that for n = 1 the minimal form factor above is 
the standard minimal form factor associated to local fields in the sinh-Gordon model computed 
in [59]. In appendix A of [27] we showed how the same expression can be derived from the 
angular quantization scheme [HH] . The branch-point twist field minimal form factors were 
also computed for the sine-Gordon model in [26] . 

3.5 Higher particle form factors 

Form factors with a number of particles higher than 2 can only be dealt with in the integrable 
case. From the consistency equations for form factors of the twist field T and generalizing the 
methods previously used for the computation of form factors of local fields in diagonal theories 
(see e.g. |59]). it is natural to make the following ansatz for the higher particle form factors: 

...,eu) = hI^'-'qI^'-\xi, . . .,Xk)llK{9^,), (3.37) 

i<j 



where 



K{e,,) = FP'\e,j)/{T), (3.38) 



and we are considering form factors of k particles of the same type (say type 1), since any other 
form factors can be obtained from these by repeated use of (|3.9p and (|3.1Up . More precisely, 

Fp^'-^'^ {ei,...,9k) = Fp^^-\ei + 2TTi{fi^ - 1),92 + 27ri(Ai2 -l),...,9k + 2m{fik - 1)), (3.39) 

with the ordering > > . . . > /^fc- is convenient to introduce the variables Xi = exp(6i/n) 
with i = I, . . . ,k. Then, exploiting the properties of the two-particle form factor discussed above, 
we find that the first two form factor consistency equations (j3.9p and (j3.10p are automatically 
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7"! 1 1 

satisfied, provided that the functions " (a^i, • • • ,Xk) are symmetric in all variables Xj. In 

7"! 1 1 

other words, " (xi, . . . , rc^) can be expressed in terms of elementary symmetric polynomials 

of the variables Xj. -ff^'^ "^ are constants. In particular, = (T) and = 1. Plugging the 
ansatz ()3.37p on the kinematic residue equations (|3.1ip and ()3.12p produces a set of recursive 
equations for the constants H^^^'"^ and the functions Q^'^ '^(xi, . . . ,Xfc) which may be solved 
for specific models. 

So far, the only model for which all solutions to these equations are known is the Ising field 
theory. The solutions were obtained in [28]. Since we are dealing the free Fermion case, it is 
natural to expect that the form factors of the twist field would admit closed expressions in terms 
of Pfaffians, as for the order and disorder fields of the Ising theory. This is indeed the case, and 
it is easy to show that 

F^^''-\ei,...,9k) = {T)Fi{K), (3.40) 

where Pf is the Pfaffian, which is mainly characterised by the property that Pf(i^)^ = det(i^) 
and K is an anti-symmetric k x k matrix, with entries 

Kij = KiOij). (3.41) 

Notice that, for the Ising model only A;-particle form factors with k even are non-vanishing. The 
function K{x) has properties 

K{9) = -K{-e), (3.42) 
K{9)\n=i = 0, (3.43) 
{K{e + is))* = -K{e-is), 6,s€R, (3.44) 

where "*" indicates complex conjugation. The Pfaffian expression is nothing else than the 
application of Wick's theorem on the particles in the asymptotic states (specialised to all particles 
being on copy 1), a contraction of two particles being K{6i2). 

3.6 Form factor consistency checks: the A-sum rule 

The form factor program above provides form factor solutions, but it is not obvious that these 
solutions are the correct ones for the branch-point twist fields T and T. Indeed, as we men- 
tioned, the twist property and invariance under other symmetries alone are not enough to fully 
fix the field; we need to specify its scaling dimension to be the minimal one. It is generally con- 
jectured that the "minimal" solution corresponds to a primary field (hence with minimal scaling 
dimension in a given Virasoro module), but it is always good to verify that the solutions found 
agree with the expected scaling dimension obtained in the ultraviolet limit. As is well-known, 
the form factor program provides a way of carrying out this verification by allowing us to com- 
pute the correlation functions of various fields of an integrable QFT. In the ultraviolet limit, it 
is possible to relate a particular correlation function to the holomorphic conformal dimension A 
(with A + A = d, the scaling dimension, and A — A = s, the spin) of a primary field by means 
of the so-called A-sum rule: 

A^ = A^ = r (e(r)t(O)) dr (3.45) 

(where the integration is on a space-like ray), originally proposed by G. Delfino, P. Simonetti 
and J.L. Cardy in [6l]- Here, where G is the local operator corresponding to the trace of the 
stress-energy tensor. The first equality, naturally expected from CFT, holds from the A-sum rule 
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thanks to the fact that G commute with T and that 0^ = G. In the cases of the branch-point 
twist fields, which are spinless, the holomorphic conformal dimension is related to the scaling 
dimension by dn = 2A-^, where dn is expected to be ()2.8p . 

Introducing the complete sum over quantum states ()3.5p and carrying out the r-integration, 
the expression above can be rewritten as 



xF®l^i-^'^(^i,...,0fc) (Fl\^^-^^{ei,...,ek)Y , (3.46) 



k ■ ■ ■ ,^ftj k 

where the sum in /ij with i = 1, . . . ,k is a sum over particle types in the theory under consid- 
eration. The sum starts at A; = 1 since we are considering "connected" correlation functions, 
that is, the A; = contribution has been subtracted. The sum above, can only be carried out in 
particularly simple cases. For most models, one must be content with evaluating just the first 
few contributions to the sum. Fortunately, the many studies carried out in the last years provide 
strong evidence that the sum above is convergent and that in fact, the first few terms provide the 
main contribution to the final result. Indeed, the convergence is often so good that considering 
only the contribution with k = 2 already provides very precise results (see e.g. [M])- Expecting 
a similar behaviour also in our case, we will approximate the sum above by the two-particle 
contribution, that is 

^ 7 7 de,de2F^^'\eu)FP'\eur 

2{T) J J 2(27r)2m2 (cosh ^1 + cosh 02)^' 

— oo — oo 

The factor of n is a consequence of summing over all particle types and using ()3.19p . In addition, 
the only non-vanishing contribution comes from form factors involving only one particle type, 
since we are considering n non-interacting copies of the model. This implies that 

F^^'^ie) = V i^j. (3.48) 

Changing variables to 8 = 6i — 02 and 9' = 9i + 62 we obtain, 

7 ,aFT{o)Fr'\er 



A' « T^-^n:^ / d9^ (3.49) 

327r2m2(T) J cosh2(0/2) ^ ^ 

—00 

In [27] we evaluated the integral above both for the Ising and sinh-Gordon models, whereas 
in [26] we performed a similar computation for the sine-Gordon model. We briefly summarise 
below our results for the Ising and sinh-Gordon theories. 

3.6.1 The A-sum rule for the Ising and sinh-Gordon models 

For the Ising model the only non- vanishing form factor of the trace of the stress-energy tensor is 
the two-particle form factor. Hence, the two-particle approximation (j3.47p becomes exact. The 
two-particle form factors are given by 

^r'\e) = . ^'!^i^.^f^L.^ ^inh ( l), F?^\9) = -2..m^sinh ( ^ , (3.50) 



^sinh(^)sinh(^) \2n 
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and therefore 



A 



T 



1 

16^ 



^°^(^)^^"Mll)^inh(f) 
sinh(^)sinh(^)cosh2(|) 



d9. 



(3.51) 



The above integral can be computed analytically for n even by shifting t by 2Trni and noticing 
that the integral changes by a sign, so that 



2A 



r 



2m rj 



1=1 



1 
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(3.52) 



where rj are the residues of the poles of the integrand at t = iTr{2j — 1), with j = 1, . . . ,n. 
These residues can be easily computed (see [27] for the details) and the expected result ()2.8p is 
reproduced, with c = 1/2. 

For the sinh-Gordon model, the relevant 2-particle form factors are given by 



Fr 



Till 



(0) 



(T)sinh(^) 



2nsinh(^) sinh 



2n 



r\ii 



(9) 



. (3.53) 



n=l 



The form factors of were computed in [59]. Since G is a local operator, its minimal form factor 
is given by (j3.36p with n = 1. The tables below show the result of carrying out the integral 
()3.49p numerically for various values of n and B. Next to each value of n in brackets we show 
for reference the expected value of A^, as predicted by the CFT formula (j2.8p (with, again, 

A^ = dn/2). 





n = 2 (0.0625) 


n = 3 (0.1111) 


n = 4 (0.1563) 


n = 5 (0.2) 


B = 0.02 


0.0620 


0.1114 


0.1567 


0.2007 


B = 0.2 


0.0636 


0.1135 


0.1599 


0.2048 


B = 0.4 


0.0636 


0.1148 


0.1620 


0.2074 


B = 0.6 


0.0643 


0.1155 


0.1631 


0.2088 


B = 0.8 


0.0644 


0.1158 


0.1636 


0.2096 


B = 1 


0.0644 


0.1159 


0.1637 


0.2098 





n = 6 (0.2431) 


n = 7 (0.2857) 


n = 8 (0.3281) 


n = 9 (0.3704) 


n = 10 (0.4125) 


B = 0.02 


0.2436 


0.2864 


0.3289 


0.3712 


0.4135 


B = 0.2 


0.2488 


0.2925 


0.3360 


0.3793 


0.4225 


B = 0.4 


0.2522 


0.2966 


0.3407 


0.3846 


0.4284 


B = 0.6 


0.2540 


0.2988 


0.3433 


0.3876 


0.4317 


B = 0.8 


0.2550 


0.2999 


0.3446 


0.3890 


0.4334 


B = l 


0.2552 


0.3002 


0.3449 


0.3895 


0.4339 



The figures obtained are extremely close to their expected value for all choices of B and n. 



4 Bulk and boundary entropy from twist field form factors 

Formulas (j2.12p and (j2.20p express the entanglement entropy in the bulk and boundary cases, 
respectively, in terms of the two-point and one-point functions of branch-point twist fields. The 
form factor expansions (j3.6p and (j3.7p then express these correlation functions as infinite series 
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obtained from form factors (j3.ip . These infinite series provide exact expressions, but also, as we 
explained in subsection 13. H the truncated, partial series give efficient large-distance expansions. 
Our main result, (jl.l2p . is an exact large-distance correction for the bulk entanglement entropy, 
and is obtained by keeping, in the series (13. 6p . only the two-particle contributions. On the 
other hand, our analysis of the boundary entanglement entropy, in the Ising model, required the 
consideration of all terms in ()3.7p . We will now explain how to go from form-factor expansions 
to entanglement entropy. 



4.1 Bulk entanglement entropy 

From expression ()3.6p and from (|2.12p . we find for the bulk entanglement entropy: 
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where C/>^°dei ^^^g defined in (|2:23|) . and 
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(4.2) 

The expression for e/c could be re-written in various ways, by performing the sum in one of 
particle indices fj,^ using invariance under shift of copy numbers, and/or by using the formula 
(j3.39p in integrable models. Also, again in integrable models, the range of rapidity integration 
can easily be extended to — oo < 6j < oo, j = 1,2, ... ,k (i.e. without ordering of rapidities), 
putting a factor l/k\. Indeed, the exchange of two rapidities in the form factors just corresponds 
to a multiplication by a two-particle 5-matrix, according to (|3.9p . By unitarity, \S^_^^^{9)\'^ = 1, 
so that the integrands in (13. 6p are the same under exchange of rapidities. Another way of 
understanding this is that with different ordering of rapidities, we simply have a different basis 
of states: instead of in-states, we have states where in- and out-configurations are mixed. We 
may also sum over these states in order to get the form- factor expansion, or we may average over 
many such bases, integrating in the full rapidity range. Such mixed bases exist as asymptotic- 
state bases in integrable models thanks to the independence of the scattering amplitudes upon 
impact parameters (see e.g. [33] )■ In non-integrable models, mixed bases are not expected to 
exist. However, thanks to (j3.13p . the latter argument can still be used, for the two-particle 
contribution, to extend the integration region to the whole rapidity range, without ordering. 
This will be useful in the next paragraph, where we study in more details the two-particle 
contribution. 



4.1.1 Next-to-leading order IR correction to the bulk entanglement entropy 



From the expression above (j4.ip . it appears that the first correction to the IR behaviour of the 
entanglement entropy would come from the one-particle form factor contribution 



ei(rm) 



lim — 

n^i an 



n 



2^ 



-rm cosh d 



(4.3) 



For all theories studied in }27l [26] this correction vanishes because their internal symmetries 
are such that only form factors of even particle numbers are non- vanishing. However, it is 
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possible to argue that even when the one-particle form factor is not zero, its contribution to the 
entropy should be. First, the one-particle form factor is always independent of the rapidity due 
to relativistic invariance and the spinlessness of the twist field. Hence, it only gives a simple 
n-dependent factor. Second, the form factor is zero at n = 1, as in that case the twist field is 
simply the identity. If we assume that it goes to zero like n — 1 (or in fact, faster than \/n — 1), 
then we can conclude that the derivative of its square at n = 1 must also vanish. Certainly, 
however, it would be desirable to evaluate exactly one-particle form factors in models where 
they do not vanish by symmetry argument, in order to verify this. 

We will now look at the next correction (the first non-trivial one), namely the correction 
associated to the two-particle form factor contribution 
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In the first line, we extended the rapidity integration range so that there is no ordering of 
rapidities (dividing by 2!) as described above, we summed over one particle index using invariance 
under shift of copy numbers, and we used formula ()3.20p . Notice that these operations are valid 
both at and out of integrability. In the second line we just carried out one of the integrals, where 
Kq{x) is a modified Bessel function of the second kind, and we define: 
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For simplicity, we assumed that there is only one particle type, so that both particles involved 
in this two-particle contribution have the same mass. As we will see, this turns out to be enough 
in order to obtain the general result. 

What follows will be devoted to show that, if we denote by f{9, n) the analytic continuation 
of f{9,n) to non-integer values of n, then: 
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Recall that f{9,l) = and that f{9,n) is only defined for integer values of n. It is also easy 
to see, from the general solution for the two-particle form factors p.28p or from the general 
expectation that the form factors vanish like n — 1 as n — > 1, that the term |i^2^'^^(^)p will not 
contribute to the derivative at n = 1. 

In order to obtain the entropy, we should now analytically continue the expression inside 
the derivative in ()4.4p . as function of rm and n, from n = 1, 2, 3, . . . to n £ [l,oo), compute 
the derivative with respect to n and evaluate the result at n = 1. The analytic continuation 
is of course not unique. The motivation for the choice of analytic continuation that was used 
in [271 [26] was based on Carlson's theorem, and on an expectation about the large-n behaviour 
of partition functions on multi-sheeted Riemann surfaces (or spaces with conical singularities). 
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Essentially, we want the analytic continuation of f{9,n), for any fixed 6, to be well-behaved 
enough at large n. In the integrable cases studied with precision in [27] and [26], we found that 
there was an analytic continuation such that the divergence at large n is less than exponential for 
any Re(n) > 0. This is a unique analytic continuation with this property by Carlson's theorem. 

Once we admit this requirement on the analytic continuation, the evaluation of the derivative 
at n = 1 can be done exactly. The observation at the basis of this evaluation is that although 
the analytic continuation of the sum in ()4.5p vanishes as n — > 1 for all 9 ^ 0, it does not vanish 
as n — > 1 when 9 = 0. That is, the function f{9, n) does not convergence uniformly in 9 as 
n ^ 1. Let us define for convenience 

/(n) :=/(0,n). (4.8) 

What we find is that its limit n — > 1 is non-zero, and positive. Note that this means that in fact, 
/(O, n) does not have a well-behaved analytic continuation from n = 1, 2, 3, . . . to n G [l,oo). 
Rather, f{n) is such an analytic continuation, but from n = 2, 3, . . . to n G [l,oo; it has the 
property that /(I) / /(0, 1) = 0. This was observed explicitly in the Ising and sinh-Gordon 
model |27j . where we found f{n) = 1/2 in the Ising model, and obtained a large-n expansion 
for f{n) in the sinh-Gordon model from which we could evaluate with some precision its value 
at n = 1. It was also observed in the sine-Gordon model ^26j, using this time a power series in 
n. See Fig. [71 

Naturally, the fact that the function f{9,n) does not converge to as n — > 1 at the point 
= is of no importance if we just want to evaluate the integral over 9 (with, as factor, a 
function that is well-behaved at 9 = 0). Indeed, the point = is of Lebesgue measure zero, 
so the value of the non- uniform limit at that point, if it is finite, does not affect the integral. 
However, we want to take the derivative with respect to n. It is possible to argue that this 
derivative -^f{9,n) at n = 1 should actually be proportional to a delta-function 5(9): it has 
non-zero support at 9 = 0. Intuitively, this is because looking at the value of f{9, n) near enough 
to = 0, we see that it should vary very fast as n — > 1, since it should go to zero for any 9^0, 
it should stay finite for ^ = 0, and it should be continuous in 9 for all n > 1. This very fast 
variation becomes an infinite variation as n — > 1 for 9 infinitesimally near to 0, which produces a 
S{9). Note that the resulting derivative at n = 1 should also be positive, since this fast variation 
goes from /(O) > to as n decreases to 1. This gives a negative correction to the saturation 
of the entanglement entropy, which is indeed expected on physical grounds. 

But why would the function f{9,n), as function of 9, not converge uniformly as n — > 1? 
In order to answer this, we have to find a way of analytically continuing in n the summation 
term in (j4.5p . There are a great many ways of doing this, and in all cases, what we find is 
that the cause for the non-uniform convergence is the collision of kinematic poles, at 9 = iir 
and 9 = 2i7rn — in, as n ^ 1, see Fig. [HI For instance, using a Poisson resummation formula 
to transform the sum over j into an integral (2^, one sees that the kinematic poles pinch this 
integral as n — > 1. Another way is to perform the summation via a contour integral with a 
cotangent kernel whose residues reproduce the terms that are to be summed. There, one has to 
extract the poles of the form factors themselves, these extracted poles collide as n ^ 1. Both 
ways necessitate periodicity of the function f{9, n) upon adding 2i7rn to 9. The latter way was 
in fact generalised to non- integrable models in |42j . 

Yet, perhaps the simplest way of extracting the derivative is simply to extract the kinematic 
poles themselves, and then to perform the summation exactly on these poles [27] . The contri- 
bution of the kinematic singularities to the sum in the function f{9, n) (j4.5p is obtained from 
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Figure 7: (a) shows 4 functions: the points are the function n/(0,n) for integer values of n in 
the interval [2, 30] both for the Ising (black squares) and sinh-Gordon (black circles) models, 
evaluated numerically. The solid green line gives the analytic continuation nf{n) for real values 
of n in the interval [0,30] for the Ising model, that is the function n/2. Notice that /(I) = 1/2 
whereas /(0, 1) = 0. Finally, the solid red line gives the function nf{oo) for the sinh-Gordon 
model, that is a straight line passing through the origin which describes the asymptotic behaviour 
of the function n/(n) for n large. In the sinh-Gordon case, all functions have been computed 
for B = 0.5. (b) is a magnification of the lower left corner of the sinh-Gordon part of (a) 
which demonstrates (14. lip , namely, as for the Ising model, /(I) = 1/2 and /(0, 1) = 0. For the 
sinh-Gordon model, the analytic continuation is n/(n) = cion + Yl^i al-i ; for some coefficients 
ao, which depend of the coupling B [2?]. Such analytic continuation was also studied 
for the sine-Gordon model in [26] and found to display a power series structure, with powers 
depending on the values of the coupling constant, very different from the one identified for the 
sinh-Gordon model. 

T"!!! / '7"lll\* 

the singular behavior in j of the summand s{9,j) = F2 (^ + 2iTij) [F2 ) — 27rij): 



i fI^^^ {29 + 2TTin - iir) i fJ^^^ (29 + iir) 

s{9,]) ~ — ^ ^ ^ + c.c. 

— 27rij + 27rin — in 9 — 27rij + itt 

(where c.c. means complex conjugate, taken for real 9 and real j). It is a simple matter to 
perform on this expression the sum "^^Z^, giving 

+ C.C. 

where ^{z) = dlogT{z)/dz is the derivative of the logarithm of Euler's Gamma function. This 
has no poles at 9 = 0, as the kinematic poles of the form factors involved that would produce 
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Figure 8: The collision of kinematic singularities when n ^ 1. 

a pole at = actually cancel out. The poles that are nearest to Re(0) = as n — > 1 are at 
9 = ibz7r(n — 1), again coming from the form factors involved. The residues to first order in n — 1 
give: 

with 

/(I) = I ■ (4.10) 

This has simple poles at 9 = ibz7r(n— 1) with residues that vanish at n = 1, gives lim^^i /(O, n) = 
/(I), and vanishes like (n — 1)'^ as n ^ 1 for 9^0. The limit n — > 1 of the derivative with 
respect to distribution on 0, is easily evaluated: 

|-/>,n)) ^= TT^f (1)6(9) . (4.11) 

The full form of f(9,n) was obtained in appendix C of ^27j. For the free case our result is in 
agreement with the n — > 1 limit evaluated previously in |65j . 

Inserting this inside (j4.ip we find that the entanglement entropy S'yi(rm) for A an interval 
of length r is: 

S\^^^(rm) = log(em) + U"''"^''^ - lKo(2rm) + O (e-^^™) , (4.12) 
3 8 

with [/model defined in (|2:23]1 . This result can be generalised to theories with more than one 
particle and/or bound states [26]: the bound states poles never collide, hence do not provide 
additional contributions; and the kinematic poles only occur for particle- antiparticle form fac- 
tors, so that we have a contribution for all particle types of the theory. The result was also 
generalised to non-integrable models: the only properties needed to derive the result were those 
seen in subsection 13.31 to hold out of integrability. Hence, we obtain the more general result 
()1.12p quoted in the introduction. 

In general, the constant [/™°dei extremely hard to evaluate, even in integrable models. 
However, for models of free particles like the Ising model, it can be evaluated by standard 
methods of angular quantisation [66\ [67] . This was done in [27] for the Ising model (by simply 
relating the branch-point twist fields to ordinary C/(l)-twist fields for which the calculation had 
already been done), where we obtained the value (|1.10p . 

A final note is in order concerning the function /(n). Above we concentrated on its behaviour 
as n — > 1, as this is what matters for the evaluation of the entanglement entropy. However, the 
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behaviour of this function as re — > oo is also somewhat interesting. Indeed, in all models studied 
where the perturbation from the UV fixed point was strictly relevant (Ising, sinh-Gordon and 
sine-Gordon models), the behaviour at large n was observed to be exactly linear. However, 
for the sine-Gordon model at a value of the coupling that makes the perturbation marginally 
relevant, we observed a behaviour like nlogn [26j. Interestingly, then, it seems that the large-re 
behaviour encodes information about the short-distance theory: the way the model behaves near 
to the critical point. 



4.1.2 Higher order IR corrections to the bulk entanglement entropy 



In order to evaluate higher order corrections to the entropy, namely 63 (mr), 64 (mr) ... we need 
first of all to know the form factors for more than two particles. As indicated in previous 
sections, such solutions are only known for the Ising model (see (I3.40p ) for which all higher 
particle corrections to the entropy were evaluated in ^28j. Once more, the problem of finding the 
right analytic continuation in re plays a fundamental role. Solving this problem has been indeed, 
even for the Ising model, extremely challenging (see [28] for the details of this). Although we 
have not investigated interacting models yet, we suspect that higher order corrections may also 
be, at least to some extent, of a universal nature, as the pole structure of the form factors plays 
again a key role. For the Ising model we have found that. 
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and e2k+iirin) = 0, where 9 
recover the result (14. 7p . 



Yliti ^iid Oij = 6i + 6j. Notice that for /c = 1 we exactly 



4.2 Boundary entanglement entropy 

In order to make use of the form factor expansion (j3.7p in the entanglement entropy (j2.20p . we 
need to describe in more detail the boundary state introduced in section 12.31 As mentioned 
there, the state \B) is just a tensor product of boundary states in the individual copies. In 
integrable models, these have an explicit expression as the famous boundary state introduced 
by Ghoshal and Zamolodchikov [29]. In the case where no boundary bound state can form, we 
have 

1^) = (^i^ E JZ^ij- ^) W j |0). (4.14) 

The function R{9) is the boundary reflection matrix of the integrable QFT (assuming still it has 
just one particle) and Z^{6) are the Faddeev-Zamolodchikov operators, which provide a gener- 
alization of the creation-annihilation operators for integrable QFTs with non-trivial interactions 
[311 168j . Their main properties are 

Z^^iOi) ■ ■ ■ Z^^iOkM = 1^1, . . . for 01 > . . . > 

z^^{e,)z,^{e2) = s,^,^{9i - 02)^M.(^2)^M,(^i). 
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The tensor-product form of the boundary state indicates that particles hving in different copies 
of the theory do not interact through the presence of the boundary. 

Using the boundary operator defined above, we can evaluate all boundary-state overlaps 
occurring in (13.70 and obtain a large-r expansion for the entanglement entropy from p.20p in 
terms of form factors, as for the bulk case. In fact, we could see this expansion simply as coming 
from expanding in (I2.20p the exponential defining the boundary state in (I4.14p . Either ways, 
this gives, with the distance between x and the boundary being r. 
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and so on. Therefore, we can write the boundary entanglement entropy as 
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(4.17) 
(4.18) 
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Here, we choose the universal model-dependent constant [/™°doi ^.j^g same as in the bulk 

expansion ()4.ip . This, as we explained in the introduction and in subsection 12.41 completely 
fixes the UV behaviour of the entropy in the boundary cases, and in particular fixes the relation 
between the non-universal short-distance cutoff e and the correlation length. 



4.2.1 Next-to-leading order IR correction to the boundary entanglement entropy 



The series (|4.19p is of immediate use to provide a large-distance expansion of the boundary 
entanglement entropy. In particular, we may easily obtain the general form of the first correction 
to saturation at large distances. From (j4.18p and (j4.20p we see that the next-to-leading order 
correction for r large is given by 
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This correction is much more model-dependent than the similar correction in the bulk case, both 
through the boundary reflection matrix and through the minimal form factors (which come from 
the scattering matrix); there is in fact little to say concerning its general features. 

In order to study the form of S2{2mr) in more detail and investigate higher order corrections, 
we now resort to the particular example of the Ising model. 

4.2.2 Boundary entanglement entropy of the Ising model 

The expansion (14.190 is not only an efficient large-distance expansion, but also an exact rep- 
resentation. Hence, it can be used to analyse the short-distance behaviour of the boundary 
entanglement entropy as well, by resumming all the terms of the series. In what follows, we will 
reproduce the main arguments and results of such a study of in the case of the Ising model, 
done in [28]. 

Let us recall the types of integrable boundary conditions that have been found for the Ising 
model, whose two-particle S'-matrix is simply —1. A family that was studied in much detail in 
[29] is that corresponding to the presence of a magnetic field that couples to the Ising spin field 
on the boundary. The spin field is the order parameter, hence we are looking at the scaling 
limit of the Ising spin chain in a transverse magnetic field whose magnitude is slightly below 
its critical value (ordered regime), and with a parallel magnetic field on the boundary. The 
corresponding boundary reflection matrix is given by 



which includes, for special values of the parameter k, the following physically different types of 
integrable boundary conditions, 

• Free boundary condition: k = 1, 

• Fixed boundary condition: k = —oo, 

• Magnetic boundary conditions (interpolating between the previous two): k = 1 — 
where h is a boundary magnetic field < h < oo. The free boundary condition would 
then correspond to /i = whereas the fixed boundary condition is equivalent to having a 
infinitely large magnetic field fixed at the boundary. 

Boundary corrections to the expectation values of the energy and disorder field in the Ising 
theory were computed using this reflection matrix in [69j . 

In the cases where k > 0, the reflection matrix has a pole on the imaginary line on the 
physical sheet, < Im{9) < in/l. This implies that the boundary state expression ()4.14p is 
not correct. A modified expression exists [29], but for simplicity, we will not analyse this case 
here. Hence, throughout we will consider k < 0. Note that the case k = does not require 
modifications, since the residue of the i?-matrix vanishes at this point. At k = 0, the bound 
state becomes weakly bound, and propagates far into the bulk 

The cases k > —1, i.e. h < he = are also somewhat special. In these cases, the 

R-matrix still has a pole on the imaginary 9 line, although not on the physical strip when k < 0. 
As noted in [29] . the case k = —1 corresponds to a "critical" value of the magnetic field, he, at 
which the reflection matrix happens to have a third order zero at = 0. 

We note also that the R matrix at k = is just equal to the negative of the fixed-boundary 
condition R matrix, k = —oo, and that 
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In order to make it clear that we are dealing with the Ising model, we will explicitly write 
the K dependence in the boundary corrections, S2fe(t) ^ S2k{t-,i^)- 



Reproducing the main arguments of [28], we will analyse the corrections S2k{t,>^) order by 

order. This will provide, by resummation of these corrections, an exact evaluation of the constant 
V{k) appearing in 

As was mentioned in subsection 12.41 we expect the constant V{k) to depend only on the 
conformal boundary condition that is reached as we bring the length r to zero (i.e. in the limit 
r <^ m^^). There are two conformal boundary conditions for the Ising model: zero magnetic 
field, i.e. free boundary condition, and infinite magnetic field, i.e. fixed boundary condition 
(although, to be precise, the latter one has two possibilities, for the two orientations of the 
boundary magnetic field). The first one is a UV boundary condition, and the second an IR 
boundary condition, from the point of view of the boundary renormalisation group flow. Hence, 
for any fixed k, we expect that in the UV limit mr 0, we obtain a situation described by the 
UV, free boundary condition. On the other hand, if the system is already at a fixed boundary 
condition n = — oo, we expect that as mr ^ we obtain a situation described by the IR, fixed 
boundary condition. Hence, V{k) should be the same for any finite k, but should take another 
value for k = —oo. 

We may wish to obtain the boundary entanglement entropy in the situation where the bulk 
is conformal, but the boundary itself is not conformal, with a magnetic boundary condition, that 
depends on the magnetic field. In this case, we have to send the bulk correlation length to 
infinity while keeping both r and the magnetic, boundary correlation length r] finite, and of the 
same order. The result will be a universal function of the ratio r/rj. This is the limit where 
is greater that all other scales, so that we must take simultaneously mr — > 0, k — > — oo in an 
appropriate way. The way this limit is taken can be understood through CFT considerations, 
but the four-particle calculation below will make it clear that it is the product Kmr that we need 
to keep fixed. In terms of the magnetic field, this means that we need to keep the product h'^r 
fixed (both ways are equivalent in the limit k — > — oo). That is, the boundary correlation length 
r] is of the order of h~^. We will not investigate this situation much further, although we will 
discuss it below in relation with the four-particle contribution. But we note that after taking 
this particular limit, if we then make the boundary correlation length rj very large or very small 
compared to r, we obtain a free/fixed conformal boundary condition. From this point of view, 
it is clear that it is indeed the point k = — oo that divides the regions in k where the UV/IR 
conformal boundary condition is reached as we take mr — > keeping k fixed. 

4.2.3 Two-particle correction to the entropy in the Ising model 

In terms of the i2-matrix ()4.22p and the minimal form factor ()3.34p , the correction ()4.2ip becomes 



We see that the correction S2{t, k) is finite for all values of t, including t = (zero distance). In 
particular, this correction is not expected to contribute to the logarithmic term in the small-mr 
behaviour of the boundary entanglement entropy, (jl.lip . At t = it is possible to evaluate the 
integral above explicitly: 





(4.24) 



, , , , 1 TT TT 

C2{k) := ^2(0,/^) = ^- g 



K (tt + 2 arcsin(K)) 



(4.25) 



4ka/1 + «^ 



and in particular 



C2(-l) 



10 - 37r 



and 



C2(0) 



7r-2 



(4.26) 



8 



8 



37 



The number c(k) will contribute to the constant V{Hi) of (jl.lip . 



4.2.4 Four-particle correction to the entropy in the Ising model 

It is interesting to analyze in some detail the next correction to the boundary entanglement 
entropy for the Ising model. For boundary theories, this is the first correction for which the 
issue of finding the right analytic continuation in n plays again a crucial role. As it turns out, 
that crucial role extends to all other higher order corrections, which were explicitly found for 
the Ising model [58]. Also, the contribution S4{t,K) gives us a lot of insight as to how the 
short-distance behaviour in (jl.lip can be recovered; in particular, how the two values of V{k) 
can come out, and how they are connected by the massless, non-conformal magnetic boundary 
situation explained in paragraph 14.2.21 above. 

The four-particle boundary correction is given by 
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Here we have used (|3.39p (in the case of the Ising model, there is no need for the ordering of the 
particle indices), as well as the notation 6^ = 9 + 27rij. Employing (I3.40p we find 
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where the function K{x) was introduced in (|3.38p . and we use the "hat" notation introduced after 
equation ()4.13p . It is simple to show that the v?K{26i)K{262) term will give no contribution 
to the derivative at n = 1 so that only the terms in the sum remain. These terms will give a 
contribution, since, employing ()3.43p - (l3.44p . they can actually be rewritten as 
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where the symbol = means equality up to integration in 6i and 62 inside the integral ()4.27p . 
We can then use the result (j4.1ip to show that the sum (j4.30p is proportional to S{9i2)- The 
four-particle correction to the saturation value of the Ising entanglement entropy is therefore 
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The four-particle contribution is divergent as t — > 0. This points to the fact that S4(t, k) will 
contribute to the logarithmic term in the short-distance behaviour in (jl.lip . Technically, the 
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reason for this is that the integrand of (j4.3ip is a function that tends to the value —1 as 9 ^ oo 
when t = 0. Therefore the integral at t = is divergent. In order to find the precise behaviour 
of the correction as t approaches 0, one can rewrite the integral above as: 
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The behaviour of the Bessel function as t goes to zero is well-known, 

i^o(2t) = -7-log(t) + 0(t'logt), 



(4.31) 
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where 7 = 0.577216... is the Euler-Mascheroni constant. Written in this form, the integral part 
is now a finite constant at t = 0, and we may define 
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with in particular 
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In view of the short-distance behaviour of the entanglement entropy (jl.lip . we expect that the 
coefficients of the logarithmic divergencies at small rm will add up to the finite number 1/12 
when all corrections are considered. This has been proven numerically in [28|. Also, the constant 
C4(k), like C2{k) above, is a part of the constant V{k) in (ll.llj) : again, in principle one should 
add up all such constants, for all corrections, in order to obtain V{k). 

Recall that V{k) takes only two possible values, one for k > —00 and one for k = —00. This 
implies that V{—oo) / limK^-00 ^(k)- It turns out that this inequality is true for C4(k), and 
should hold as well as for infinitely many constants C2ii{k) [2H]. Technically, we observe that 
although the integral in (I4.3ip has the same value for k = — 00 as for k = for any t > 0, 
we have lim^ — ^—00 

C4(k) = 3/8, different from 04(0) = 1/8. The explanation is that the hmit 
t ^ of the integral in ()4.3ip as a function of k is not uniform. For all values of t > we have 
Si{t, —00) = S4(t, 0), and there is a maximum for k S (—00, 1) at a unique value k = kq. But 
as t becomes smaller, the position of this maximum shifts towards more negative values, until 
it reaches —00 at t = 0. There, if we take away the constant (as function of k) term log(t) 
in order to make the limit finite, the value of the maxiniuni itself readies lini/^^—oo 

C4(k). It is 

also possible to observe in the integral in (I4.33P that the symmetry between k = —00 and k = 
is broken. Indeed, if k ^ —00, the term in parenthesis can be approximated by 1 except for 
values of 9 where k + cosh^ ~ 0. These are very large values of 9, but they are not damped by 
any other factor, hence the mistake in approximating by 1 is non-negligible for any k. 

This means that the expansion (I4.35P is valid only for k > —00. For the case k = —00, 
that is, the fixed boundary condition, we have to consider the other order of the limits: first 
K —00, then t ^ 0. By the symmetry between k = —00 and k = 0, we define 
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so that ()4.35p still holds in the case of a fixed boundary condition, k = — oo. 

Recall our discussion in paragraph 14.2.21 about the particular massless limit whereby the 
boundary is still magnetic, not conformal. In this perspective, it is instructive to obtain a more 
general small-t expansion, where we take simultaneously k — > — oo. Let us consider t — > with 
—nt = a fixed. We may use the change of variable s = cosh0 — 1 and write S4(t, k) as 

IQJo \k-l-sj V (s + 2)3/2 s^ij 16 \k-l-sj s + 1 

(4.37) 

The first integral as a function of k has a uniform limit as t ^ on k G [— oo, 0), so that we can 
directly take k = — oo and t = 0; this gives (2 — log2)/16. The second integral does not have a 
uniform limit, but it can be evaluated explicitly: 

4ae"2* - (a + t)r(0, 24)6^* - 8a{a + t)e'^^^+^^r{0, 2(a + t))) 

where T(z, u) is the incomplete Gamma function, v^~^e~'"dv. The small-t limit can then 
easily be taken: 

S4(t, -a/t) = ^ log(t) + + 4(a) + 0(t). (4.38) 
lb lb 

where 

4(a) = ^-iae2-r(0,2a). (4.39) 
It is easy to see that c\{a) interpolates between limK^-oo C4(/i) at a = to C4(— oo) at a = oo. 

4.2.5 Higher order corrections to the entropy in the Ising model 

Many of the subtleties observed in the previous section for the four-particle correction to the 
boundary entanglement entropy generalise to higher orders, in particular the need to find the 
analytic continuation in n of our expressions and the non-commutativity of the limits observed 
above. We do not wish to go into the details of the computations involved, which are very 
technical and cumbersome, even for the Ising model. However, let us give an indication of how 
it works. 

The main complication arises from the fact that the pole- and zero-structure of the form 
factors (in particular, the way some poles and zeroes cancel each other) changes substantially 
as soon as n is allowed to take non-integer values. As a consequence, the phenomenon observed 
in [2Tj, and recalled above, of non-uniform convergence of form factors as n ^ 1 is generalised 
to a non-uniform convergence as n — > A;' for all positive integers k' < k for the 2/c-particle form 
factor. In order to obtain the correct analytic continuation, the principle is that it is obtained 
from the analytic function that describes form factor contributions at values of n large enough. 
This amounts to evaluating first the contribution we would obtain analytically continuing in 
n around some integer and then adding the residues of all the extra poles that are crossed by 
the integration contours when bringing n from infinity. This is naturally not the only analytic 
continuation that is possible and therefore it becomes quite crucial to find ways of checking it 
for consistency. The exacy UV behaviour shown below provides good support, see [28j for a 
more extensive discussion. 
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The general formula for all higher-particle corrections found in |28j takes the form: 
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and similarly 
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both formulae hold for k = 1, 2, 3, . . .. All integrals involved are absolutely convergent. Notice 
that for /c = 1 in (|4.4U|) we recover the result (j4.3ip as it should be. Setting A; = in (|4.41|) does 
not give (I4.24p . since /c = is out of the range of applicability of these formulae; the two-particle 
case S2{t, k) is a special case. 

From these expressions it is also possible to extract the divergent part as t ^ 0, as done 
for the four-particle contribution in the previous paragraph. This is achieved by simply setting 
t = in all expressions, after subtracting 1 from the product of i?-matrices. We then obtain 
S4kit,K) ~ C4fc(K) + e2k{t) and Sik+2{t,K) ~ Cik+2 as t ^ 0, for all k = 1,2,3, .. ., where e2k{t) 
are the expansion coefficients ()4.13p in the expression for the bulk entanglement entropy (see 
(j4.ip . and where 
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and similarly 
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both formulae for A; = 1, 2, 3, . . .. The constant V{k) that characterises the UV behaviour of the 
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boundary entanglement entropy as shown in (jl.lip can now be obtained as follows: 
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hence 



V{K) = Y,C2k{K). (4.45) 

fc=i 

Naturally, in order for ()4.45p to be a correct representation of V{k), the infinite sum over k 
should give a finite result. This is quite subtle, as form factor expansions are expected to 
provide convergent series expansion for finite distances, but not necessarily at zero distance. In 
[28j, exact evaluations of the first few coefficients C2fe(K) for k = 0, — 1,— oo, and extrapolation 
to higher k, gave strong indications that the series is indeed convergent. More precisely we have 
managed to obtain closed formulae for C2fc(0) and C2fc(— cxd): 

1 IT / 2k — 1 \ ^ 1 

'''^^^ = 8k{2k - ly ^^'^^^(o) = ^ k-1 ) -M2kTTy ^^-^^^ 

and C4fc(— oo) = C4fc(0) and C4k+2{—^) = — C4fe+2(0). This then gives 

2 °° 

V{0) = ^ + C4fc(0) = log V2, (4.47) 
fc=i 

and 

2 — vr °° 

y(-oo) = ^— + ^(C4fc(0) - C4fc+2(0)) = 0. (4.48) 

k=l 

We have also computed the values of C2fc(— 1) for /c up to 4 and found that V{—1) ~ 0.321966.... 
This value is compatible with log-v/2 = 0.346574.... Therefore, as expected, in the UV limit the 
entropy only depends on whether the magnetic field h is finite (free boundary conditions) or 
infinite (fixed boundary conditions). We can then summarise our result as 

[0 [K = — oo). 

4.2.6 Connection with the boundary entropy 

It is natural to interpret V{k) as a "boundary entanglement": the contribution of the boundary 
to the entanglement between the region A and the rest. Naturally, for fixed boundary condition, 
there should be no contribution at all, since the boundary does not experience quantum fluctu- 
ations. Our result ()4.49p shows that we have chosen the correct large-distance normalisation to 
have V{—oo) = 0. On the other hand, for free boundary conditions, the boundary fluctuates and 
should participate to the entanglement. This is in agreement with V{k > — oo) = log \/2 > 0. 
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In fact, it was shown in |28j that we may connect V{n) to the boundary entropy s, a quantity 
that essentially counts the number of degrees of freedom pertaining to a boundary. This quantity 
is simply given by s = log 5 where g is the boundary degeneracy introduced by Affleck and 
Ludwig jH]. In particular for a bulk CFT, they showed that g = {0\B) where |0) is the bulk 
CFT ground state, and \B) is a boundary state of norm 1 in the bulk CFT Hilbert space (in 
particular, s < 0). It is considering the boundary entropy that Friedan and Konechny [70] were 
able to provide a proof of the "(7-theorem" : that the (7-function decreases in the RG flow from 
UV to IR. For the Ising model, it turns out that s = in the free boundary case, and s = 
in the fixed boundary case. 

In order to understand the general relation between V{^k) and s, we must make use of the for- 
mula derived by Calabrese and Cardy [,22j connecting the boundary entropy to the entanglement 
entropy at criticality: 

S = S^°"'^"^^^(r)entical " ^5^"'^2r)„itical. (4.50) 

The result is independent of r, since both entanglement entropies have the same logarithmic 
r-dependence. Naturally, the universal entanglement entropy, as we emphasised in the intro- 
duction, must be defined by subtracting a non-universal constant, the logarithm of the corre- 
lation length. Hence, the universal entanglement entropy is only defined up to addition of a 
r-independent and boundary-independent constant. We provided in the introduction precise 
ways of fixing this constant in the bulk and boundary cases, by looking at asymptotics, with 
a massive bulk. However, formula ()4.50p holds if and only if both entanglement entropies are 
evaluated in the same cut-off scheme. For instance, both should be evaluated on the same lattice 
at infinite correlation length, and taking the distance r to be the same number of lattice sites 
(this number tending to infinity). Formula (|4.50p is not of immediate use in the context of QFT, 
because the requirement of a consistent cut-off scheme is a not a universal requirement. 

It turns out that the use of massive QFT solves this problem. More precisely, we can make 
sure that we have the same cut-off schemes in calculating 5^°"°*^^'^^(r)criticai and >S']4""^(2r) critical 
by connecting them to our universal definitions of the entanglement entropy in the bulk and 
boundary massive cases. Recall that our definitions involve the inherently massive large-distance 
asymptotics. More precisely, we can show that, in the Ising model, 

V{k) = 5^°'^°^"'^^(r)criticai - ^5r'^(2r)criticai + log ^2. (4.51) 

The entanglement entropies should be evaluated at the UV fixed point if k is finite, and at the 
IR fixed point if k = — oo. 
This implies that 

y(K) = s - log \/2. (4.52) 

In fact, it was argued in [28] that in general the shift log \pl should be replaced by — log C where 
is the fraction of the ground state degeneracy broken by the boundary condition for large 
enough h (at the IR point). In the Ising model, the ground state has a double degeneracy (the 
two orientations of the spin), and a large magnetic field h breaks it to 1, so that = 1/2. This 
implies 

y(K) = s-logC. (4.53) 

Such an "extra" contribution — log C to s = log g was also found in [7T] in calculations using 
Thermodynamic Bethe Ansatz techniques (although it is important to note that the Thermo- 
dynamic Bethe Ansatz used there is fundamentally different from our approach based on form 
factors). In [71] the fiow of g between critical points was studied for several families of minimal 
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Toda field theories with ground state degeneracy k and with a boundary completely breaking 
it. The factor C = was termed "symmetry factor". Interestingly, it was shown that at the 

infrared point of these models, one gets g = C. Hence, for massive models with spontaneously 
broken order-parameter symmetry, and with an order-parameter boundary perturbation, one 
should find V{k) = s — sjh > 0, where sm = logC is the infrared value of s. For the same order- 
parameter perturbation both on the bulk and boundary, one should simply find V{k) = s < 0. 
All these considerations should not depend on integrability. 

Let us now discuss briefly how we obtain (|4.51|) . following [2H]. Consider the entanglement 
entropy S^^'^'^'^'^^{xi,X2) defined in subsection 12.31 but with a slightly different normalisation 
specified below; we will denote it SAixi,X2)- We may uniquely fix the cutoff, for instance by 
requiring the conformal normalisation, SaIxi, X2) ~ | log(|a:2 — xi\/e) + o(l) as X2 ^ xi. But 
this short distance behaviour is the bulk critical entanglement entropy, so that we have chosen 
e with 5'^'^^'^(r)criticai = §log(?'/e)- Using the same object, hence with the same lattice spacing, 
we can also define the critical entanglement entropy in the boundary case, 5'^°"°'^^'^^ (r) critical) 
following subsection 12.31 Note first that for both X2 and xi far from the boundary and far from 
each other, the entanglement entropy saturates to some constant — |log(me) + U thanks to 
the presence of the mass. This saturation is a sum of the contributions of the two boundary 
points at xi and X2, each one contributing — | log(me) -|- U/2. Now for X2 alone being far from 
the boundary, we get our usual boundary case, but with an extra contribution of the boundary 
point X2 at infinity. Subtracting this contribution, we find SAir) with r the distance between xi 
and the boundary: this is our usual boundary case, with possibly a different cut-off definition. 
Taking the limit rm 0, we get ^^(r) ~ 5'^°''°'^''''^(r-)criticai + o(l). 

We must now relate SA{r) to 5^°"'^'^'^^^ (r) . We know that for r large, the former saturates 
to — |log(me) -|- U/2, and the latter to — ^log(me) + U/2. The difference between U and U 
is that the former really just counts the contributions to the entanglement entropy around the 
boundary points of A when A is large, whereas the latter counts an additional entropy due to 
the ground state degeneracy. Indeed, the constant U occurs in the bulk, which should be seen as 
a large-volume limit of a periodic space where both spin-up |0, -|-) and spin-down |0, — ) ground 
states are counted. The constant U occurs when the system has a boundary, where the ground 
state is fixed if the magnetic field is large enough. Hence we have 

U = U + log2. (4.54) 

Then we find S^°'"''^'"'^(r)criticai = f log(r/e) + V{k) - log \/2, which shows (gSI]). 

There is a nice technical way of obtaining (j4.54p . using branch-point twist fields [28j. Al- 
lowing fluctuations amongst the two degenerate ground states corresponds to choosing |0) = 
((|0, -|-) -|- |0, — ))/-v/2)®"' in the n-copy model. Our branch-point twist field form factors assume 
the use of this ground state. This can be seen by recalling that particles in the Ising model 
are domain walls, separations between up- and down-spin domains. In order to have a non-zero 
matrix element of the branch-point twist field between a state with a particle on one sheet and 
another particle on another sheet, and the n-copy ground state, there must be components of 
the ground state where one sheet is on |0,+) and another sheet is on |0, — ). See figure O A 
completely up ground state |0, -|-)®", for instance, wouldn't admit such form factors. Since all 
configurations of an even number of particles distributed amongst the n sheets are allowed, we 
must use the completely symmetric ground state. Then, the constant U is the large-distance 
limit of the bulk entanglement entropy evaluated from (0, -|-|®'^T(ri)T(r2)|0, -|-)®", whereas U 
is that obtained from (0|T(ri)T(r2)|0) with the symmetric ground state. The difference can be 
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Figure 9: How a vacuum of the form • • • (X" |0, +) (8) |0, — ) (8) |0, +) (8) • • • becomes a two-particle 
state after going through a branch-point twist field. 

computed explicitly. Note that matrix elements 

(0, e'll • • • (0, e'jr(ri)r(r2)|0, ei) • • • |0, e„) 

have zero large-distance limit unless Ci = ej = e'j for all i,j (that is, all signs are the same), since 
otherwise domain walls will have to propagate between the twist fields. Hence we immediately 
find 

(0|T(ri)t(r2)|0) ~ 2'-^\0, +|^"T(ri)t(r2)|0, +)®" (4.55) 
at large distances, so that, taking derivatives with respect to n, we have 

- (|^(0|T(ri)t(r2)|0)) ^ log 2 - (^|-(0, +r"T(ri)t(r2)|0, +)®") (4.56) 

which indeed gives (|4.54|) . 

5 Conclusions and discussion 

In this review we have summarised the fundamentals and main results of a program which was 
proposed in [2^ and developed further over a series of recent publications [Ml SSI El]- This 
program allows for the computation of the entanglement entropy of a connected region A a 
massive quantum one-dimensional system with respect to the remaining part of the system. To 
this aim, it takes full advantage of integrable quantum field theory techniques, as well as analytic 
scattering theory beyond integrability, by relating the bi-partite entanglement entropy to the 
correlation functions of a particular kind of local quantum fields: branch-point twist fields. 

More precisely, we take as starting point the well-known "replica trick" and realise that, 
starting with a quantum integrable model and constructing a new model consisting of n non- 
interacting copies of the original theory, a new local QFT is obtained which naturally possesses 
permutation symmetry. Associated to certain elements of this symmetry, two twist fields T and 
T exist, whose correlation functions are directly related to the bi-partite entanglement entropy. 
In general, there are as many branch-point twist fields in a correlation function as there are 
boundary points of the region A. We considered two cases, where the whole system does not 
have boundaries (bulk two-point functions), and where the system has one boundary and the 
region A starts at this boundary (boundary one-point function), see Fig. [TJ Through the replica 
trick, the entropy is the derivative with respect to n of the correlation function evaluated at 
n = 1. 
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Since the fields T and T are local fields of the n-copy theory, and the boundary state 
admits a realization in terms of the scattering states of the bulk model, the correlation functions 
(two- and one-point functions) above can be computed by using the form factor expansion of 
massive quantum field theory. They can be expressed as a sum for different particle numbers 
over products of the form factors of the fields involved. This gives a large-distance expansion: 
computing the correlation function in the two-particle approximation, for instance, gives the 
behaviour at large distances. This expansion is in fact expected to converge rapidly, and the 
two-particle approximation is often enough to describe the correlation function up to relatively 
small distances. 

The leading behaviour of the entropy as a function of the distance was already well-known 
for very short and very large distances, both for bulk and boundary models, before our work and 
is described by the equations (jl.9p and (jl.lip . For r <^ m^^ (m being the mass of the lightest 
particle) that behaviour is fixed by the properties of the underlying CFT which describes the 
integrable QFT in the ultraviolet limit (in particular, its central charge). Thus, the entropy can 
be computed explicitly by using CFT techniques [22l[2l]. At large separations r » (in the 
infrared limit) the entropy is known to saturate to a constant value. 

One of the main results of our work has been to evaluate the first correction to the bulk 
entropy in the infrared (large-distance) limit, providing therefore a description of the behaviour 
of the entropy in the intermediate region of values of rm. This correction is obtained from the 
two-particle contributions to the form factor expansion. In fact, the two-particle approximation 
provides both the saturation value of the entropy, coming from the disconnected part of the 
correlation function (the square of the vacuum expectation value of the twist fields), and the 
exact first corrections up to 0(e~^''™), coming from the two-particle contributions. The exact 
value of the saturation in the Ising model (ll.lOp was computed in [27] and showed there to be 
in good agreement with previous numerical results |13j . 

The most surprising result of this analysis has been to establish that the leading correction 
to the entropy at large rm is in fact a more universal quantity than expected, that is, it does 
not depend on the particular scattering matrix of the model we started with, but only on the 
spectrum of masses of the particles of the original theory. This was shown to be true for all 
integrable QFTs (with or without backscattering and/or bound states) |27t and even for 
non- integrable ones [l2]. It is quite remarkable that the entropy should encode so explicitly 
crucial information about the theory both in its UV regime (the central charge) and in the IR 
regime (the number of light particles). 

We have deduced this result from general arguments and checked it explicitly for the Ising, 
sinh-Gordon models and sine-Gordon models. The mathematical reason for this "universal 
behaviour" is clear, as the result is directly related to the presence of kinematic poles in the 
two-particle form factors. Only form factors having such poles do contribute to the final result 
for the entropy, and their individual contributions turn out to be independent of the scattering 
matrix. The presence of bound state poles does not change this conclusion, neither does the loss 
of integrability. 

For the quantum Ising model, we extended the main result above and identified exact infinite- 
series formulae for the bi-partite entanglement entropy both in the presence and absence of 
boundaries |2H]. In order to obtain our formulae, we found closed expressions for all non- 
vanishing form factors of branch-point twist fields in the n-copy theory; we identified the correct 
analytic continuation in n of the contributions of these form factors to correlation functions and 
evaluated their derivatives with respect to n; and we checked both the form factor formulae and 
their analytic continuation for consistency. Since Wick's theorem applies, all form factors of 
the twist field admit expressions in terms of a Pfaffian. From this, obtaining the right analytic 
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continuation of every contribution to the form factor expansion of the twist-field boundary 
one-point function or bulk two-point function is a highly non-trivial problem. We solved this 
problem in [28j, and verified it for consistency in a very precise manner by finding an explicit 
formula for the leading logarithmic behaviour both of the two-point function of the twist field 
in the bulk and of its derivative at n = 1, and using a combination of analytical and numerical 
computations. 

The results just described for the Ising model have put us in the position to analyse another 
quantity of interest in this context, that is the contribution to the free energy that can be 
attributed exclusively to the presence of the boundary. This is essentially the boundary entropy, 
the natural logarithm of the boundary degeneracy or ^(-factor originally introduced by Affleck 
and Ludwig in [44j- In our analysis we have computed the universal quantity V{k) which is 
closely related to g (where k is related to the boundary magnetic field in the Ising model). We 
have defined V{k) in a universal, QFT way as a certain rm independent contribution to the 
boundary entanglement entropy in the UV limit. We have found an exact formula for V{k,) and 
evaluated it at k = 0, — oo to V{0) = log \/2 and V{—oo) = 0. We have also gathered strong 
numerical evidence that V{k) is in fact constant and equal to V{0) for any finite values of k. 
These two values of V{k,) would correspond to the two conformal invariant boundary conditions 
that are known for the Ising model: the free and fixed boundary conditions. The fact that V{k) 
is larger for free boundary conditions (finite magnetic field) and that V{0) — V{—oo) = log-v/2 
are properties which also hold for log g. However, it is known from Cardy and Lewellen's work 
[50] that (7fi.ee = 1 and ga^cd = 1/ V^, hence 

^(k) -logc/ = log\/2. (5.1) 

We have identified the difference between these two quantities as an IR contribution to the 
entanglement entropy coming from the ground state degeneracy in the periodic case that is 
broken in the boundary case. We have proposed a generalisation of this to more general models 
with relevant boundary perturbations. 

Although not discussed here in any detail, we have also found interesting results for the IR 
behaviour of the entanglement entropy in the boundary case. It is known that the bulk entropy 
saturates for large distances. In particular, the fact that the entropy is an increasing function 
of rm follows from the "strong subadditivity theorem" and translation invariance, as proven in 
|17tl72j. Since translation invariance is broken in the boundary theory, the entanglement entropy 
in this case is not necessarily a monotonic function of rm. Indeed, we find that, for a range of 
values of k, it has a maximum for some value of rm before reaching its asymptotic value. This 
range of values of k seems to start at k = — 1 which corresponds to a value of the magnetic field 
for which, in a sense, the boundary becomes "critical" . 

There are many open problems related to our work and in general, to the computation of 
the entanglement entropy in integrable QFT. In the case of the bulk Ising model, it is known 
that the entanglement entropy can be described via Painleve transcendents [651 ES] ■ It would be 
interesting to check the consistency of this representation with our full form factor expansion. 
In the boundary Ising case, it is known that the one-point function of the order parameter has 
a Fredholm determinant representation for any magnetic field, from which differential equations 
can be derived [69]. It would be interesting to see if similar formulae hold for the branch- 
point twist field. In the general QFT case, the most obvious problem is perhaps to extend the 
present boundary analysis to theories other than the Ising model. We believe that this should 
be possible to some extent but it is very unlikely that re-summations can be done analytically 
for interacting models. Yet an independent check of (14.530 in the more general situation would 
be useful. It would also be interesting to apply the form factor approach to the computation of 
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the entanglement entropy of disconnected regions (see e.g. [THITS]), both for bulk and boundary 
theories and also to extend the analysis of the present paper to the finite temperature situation. 

An additional result of our work has been to develop the form factor program for branch-point 
twist fields in integrable models. As a consequence of the particular exchange relations between 
these fields and the fundamental fields of the theory, the form factor consistency equations 
for branch-point twist fields are different from those associated to standard local fields. In 
particular, the crossing and kinematic residue equations are modified. Since here we have only 
been concerned with the two-particle form factors of branch-point twist fields, an interesting 
open problem remains, namely to find closed solutions to the form factor equations for arbitrary 
or at least higher particle numbers for models other than the Ising theory. As explained in 
section [331 we expect these solutions to be given in terms of elementary symmetric polynomials 
of the variables e^'/". In fact, the form factor program for more general twist fields is an aspect 
of integrable QFT that still needs much more development. 

Acknowledgments: We are grateful to J. L. Cardy for his collaboration in the paper |27j . 
upon which this review is partly based. We also thank D. Bernard for pointing out to us the 
paper [19] . 
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